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ABSTRACT 

We study the excitation and damping of tides in close binary systems, accounting for the leading 
order nonlinear corrections to linear tidal theory. These nonlinear corrections include two distinct 
physical effects: three-mode nonlinear interactions, i.e., the redistribution of energy among stellar 
modes of oscillation, and nonlinear excitation of stellar normal modes by the time- varying gravitational 
potential of the companion. This paper, the first in a series, presents the formalism for studying 
nonlinear tides and studies the nonlinear stability of the linear tidal flow. Although the formalism 
we present is applicable to binaries containing stars, planets, and/or compact objects, we focus on 
non-rotating solar type stars with stellar or planetary companions. Our primary results include 
the following: (1) The linear tidal solution almost universally used in studies of binary evolution is 
unstable over much of the parameter space in which it is employed. More specifically, resonantly 
excited internal gravity waves in solar-type stars are nonlinearly unstable to parametric resonance 
for companion masses M' > 10 — lOOM® at orbital periods P w 1 — 10 days. The nearly static 
"equilibrium" tidal distortion is, however, stable to parametric resonance except for solar binaries 
with P < 2 — 5 days. (2) For companion masses larger than a few Jupiter masses, the dynamical tide 
causes short length scale waves to grow so rapidly that they must be treated as traveling waves, rather 
than standing waves. (3) We show that the global three-wave treatment of parametric instability 
typically used in the astrophysics literature does not yield the fastest growing daughter modes or 
instability threshold in many cases. We find a form of parametric instability in which a single parent 
wave excites a very large number of daughter waves [N sa 10'^[P/ 10 days] for a solar-type star) and 
drives them as a single coherent unit with growth rates that are a factor of ~ TV faster than the 
standard three-wave parametric instability. These are local instabilities viewed through the lens of 
global analysis; the coherent global growth rate follows local rates in the regions where the shear 
is strongest. In solar-type stars, the dynamical tide is unstable to this collective version of the 
parametric instability for even sub- Jupiter companion masses with P ;^ a month. (4) Independent 
of the parametric instability, the dynamical and equilibrium tides excite a wide range of stellar p- 
modes and g-modes by nonlinear inhomogeneous forcing; this coupling appears particularly efficient 
at draining energy out of the dynamical tide and may be more important than either wave breaking 
or parametric resonance at determining the nonlinear dissipation of the dynamical tide. 
Subject headings: binaries: close - hydrodynamics - planetary systems - stars: interiors - stars: 
oscillations — waves 



1. INTRODUCTION 

In binary systems at small separation, the differential 
gravity from one body can induce a significant tide in the 
other. As long as the tidal potential is time- varying, the 
tide and orbit exchange energy and angular momentum. 
Dissipation damps away the variable component of the 
tide, leading to an evolution of both the orbit and the 
spins of the bodies. The final state of such dissipation, if 
it is allowed to continue, is a circular orbit, with the spin 
and orbital angular momentum aligned, and the bodies 
synchronized with the same side always facing each other. 

The rate at which the orbit and spins evolve depends 
on two quantities: the amplitude of the tide and how 
rapidly the energy in the tide is dissipated ("tidal fric- 
tion"). By analogy with a damped, driven oscillator, 
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energy dissipation at a rate E acting on a tide with fre- 
quency Lo and energy E induces a phase lag 5 ^ E/{u]E) 
between the time of maximum tidal acceleration and high 
tide. This phase lag is often written in terms of the tidal 
quality factor Q ~ 1/(5, where large Q implies little dissi- 
pation a nd vice v ersa. First modeled in the classic tre at- 
ment of IDarwinI (1879; see IGoldreich fc Soteil [l966l for 
applications in the Solar System) , tidal friction is funda- 
mental to our understanding of the origin, evolution, and 
fate of a wide variety of binary systems, including plan- 
ets and moons in the solar system, extrasolar planets, 
solar-type binaries, and compact object binaries. 

What is the physical origin of the 'friction' which cir- 
cularizes the orbits and synchronizes the spins of com- 
pletely fluid bodies in binary systems? Most investiga- 
tions treat the fiuid motion as a linear perturbation of the 
background state, and focus on linear damping mecha- 
nisms such as radiative diffusion or viscosity from eddies 
in convection zones. However, these applications of lin- 
ear theory fail to explain many of the observed properties 
of binary systems, as we describe below. 

This paper goes beyond the standard linear theory 
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treatment in order to understand the effects of non- 
linear interactions on tides in binary systems. Al- 
tliougli asp ects of this problem Iiave been studied prcvi- 
ously (e.g.. iPress et al.lfTOTSt iKumar fc G oodman 1996; 
iGoodman fc D icksonlll998t iBarker fc Ogilvia 2010), we 
present a comprehensive formalism for understanding 
nonlinear tides in stars, planets, and even compact ob- 
jects. The ideas and formalism we present are general, 
but for concreteness all of our examples focus on tides in 
slowly rotating solar type stars with either stellar or plan- 
etary companions; in such stars internal gravity waves 
(g-modes) can be nearly resonant with the tide. Our 
goal is ultimately to quantify the orbital and spin evo- 
lution of binaries accounting for nonlinear damping and 
excitation of tides. In this first paper, however, we focus 
on determining the conditions under which standard lin- 
ear theory is invalid because nonlinear instabilities drain 
energy out of the linear tide - we leave calculations of the 
resulting Q-values and orbital evolution to future work. 
Our study of tides is motivated by a number of obser- 
vational puzzles, including: 

1. Orbital circularization of solar binaries — Binaries 
containing two solar-type stars, or a solar-type star 
and a gas giant exoplanet, are observed to circularize 
much more rapidly than expected from even the most 
detailed theoretical stud ies (e.g., lOgilvie fc LinI 120041 : 
iMeibom fc Mathieul 120051 ). The observed circularization 
rates imply Q ~ 10^^ - 10^, a factor of 100 - 1000 times 
smaller than theoretical predictions. These observations 
test theories of tides in solar-type stars and gas giant 
planets over a range of orbital periods and amplitude of 
the tide. 

2. Orbital decay of hot Jupiters — The tidal Q inferred 
from observations of the circularization of solar-type bi- 
naries makes dire predictions for the orbital decay of hot 
Jupiters (Jupiter-size extrasolar planets with orbital pe- 
riods < 1 week). Tides raised in the star by the planet 
attempt to spin up the star, decreasing the orbital semi- 
major axis to conserve total angular momentum. If the 
tidal Q for solar-type stars is really 10^ — 10®, Dar- 
win's theory predicts that all 5 Gyr old hot Jupiters 
insid e 4 day orbits shoul d have spiraled into their host 
(e.g.. lOgilvie fc Lin|[2007[) : by contrast, there are -- 100 
such planets observed. The observed population of hot 
Jupiters could represent the tail end of a much larger 
parent popula tion, most of w hich have decayed into their 
host stars (Jackson et al.|[2009) . Alternatively, the (near 
universal) practice of calibrating Q from one observa- 
tion and applying i t to another may be incorrect (e.g., 
lOgilvie &: LinI [20071 ): it is certainly suspect without a 
deeper understanding of the underlying physics. For ex- 
ample, such a calibration would be invalid if tidal excita- 
tion and dissipation depends on orbital frequency and/or 
the amplitude of the tide (i.e., the mass of the secondary). 
In this paper, we show that this is indeed the case. 

3. Radii of hot Jupiters — Tidal heating has been in- 
voked to explain the large observed radii of transiting 
hot Jupiters compared to evol u tionary mo dels (Gu ct al.! 
[20031: lArras fc Bildste^ [20061: I Arras fc Socrates 2010) . 
Primordial eccentricity and asynchronous rotation will 
give rise to tidal heating, the duration of which is set by 
the value of Q. Steady state mec hanisms, such as ther- 
mal tides (jArras fc Socrate"il2010D , may also be at work. 



for which the tidal heating rate is also dependent on Q. 
A question which has received less attention is the depth- 
dependence of tidal heating. In order for tidal heating 
to effectively bloat radii, the heat must be deposited in 
deep regions with long thermal times. Understanding the 
depth dependence requires detailed calculations of wave 
excitation and damping. 

The observational puzzles above, and the physical mo- 
tivation that nonlinear effects can be important, led us 
to the present investigation. At first glance, it may be 
surprising that nonlinear effects are important, when ba- 
sic estimates such as the height of the tide relative to 
the stellar radius are small. However, the importance of 
nonlinear effects is more subtle since wave amplitudes can 
become large in highly localized regions. Consider, e.g., 
surface waves on the ocean which can travel great dis- 
tances relatively undamped, but are susceptible to wave 
breaking as the depth decreases near shore. 

In the next subsection we review the results from linear 
tidal theory in more detail; we also relate our nonlinear 
calculations to these standard linear theory results. In 
§[121 we summarize the structure of the remainder of this 
paper. Throughout we presume a basic understanding 
of the physics of stellar g-modes and p-modes (see, e.g., 
lUnno et al]|1989l:lChristensen-Dalsgaardll2003[ ). 

1.1. Linear Theory Predictions for Solar- Type Stars 

The tidal response of a fluid is typically decomposed 
into "equilibrium" and "dynamical" components. The 
simplest approximation for the tidal flow, the "equilib- 
rium tide," ignores fluid inertia, allowing the fluid to fol- 
low gravitational equipotentials. The equilibrium tide 
has large lengthscales, and hence is weakly damped by 
diffusive effects (such as turbulent viscosity due to con- 
vection), leading to ci rcularization times m uch longer 
than are observed (e.g. . IGoodman fc Ohlll997l) . The "dy- 
namical tide" refers to the resonant excitation of waves 
by the tidal acceleration. This, together with linear 
mechanisms of dissipation, is the key physics in most 
linear predictions of tidal dissipation. 

In this paper, we calculate the nonlinear stability 
of both the equilibrium and dynamical tides. The 
stability of the equilibrium tide in particular has 
rarel y been considered in the literature (although see, 
e.g.. Press et al. '1975', 'Pap aloizou fc Pringld [l98ll and 
;Kumar & Quataert 1998 for some discussion along these 
lines). It is, however, an a priori promising source of en- 
hanced orbital evolution given that the equilibrium tide 
contains most of the energy in linear theory. 

Solar-type stars are composed of a stably stratified core 
and a neutrally stratified convective envelope. Ignoring 
stellar rotation, tides in the convective envelope resem- 
ble the equilibrium tide. However, buoyancy in the core 
gives rise to short wavelength internal gravity waves that 
can have periods comparable to the orbital period of the 
binary; these waves have large amplitudes at the stel- 
lar center, but small amplitudes in the convection zone. 
Lastly, including both stellar rotation and fluid inertia, 
the core gravity waves are altered, and inertial waves 
restored by the Coriolis force are possible in both the 
radiative and convective zones. 

The main linear dissipation mechanisms for solar-type 
stars are radiative diffusion in the core and turbulent vis- 
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cosity due t o convection in the envelope . iTerguem et aLl 
([199^ and iGoodman fc Dicksonl (|1998[ ) carried out de- 
tailed calculations for nonrotating stars. Both studies 
show that resonant gravity waves in the core are only 
weakly damped by radiative diffusion, giving Q values far 
too large to be of interest. Goodman & Dickson sharp- 
ened the question. They showed that even if the train 
of gravity waves traveling inward from the core-envelope 
boundary - where they are excited by the tidal potential 
- were completely absorbed at the center the resulting 
dissipation would be Q ~ 10^ [P/(10 days)]®/'^ (based on 
their eq. [13]), too large by 3 orders of magnitude. Good- 
man & Dickson did not present detailed calculations of 
the damping of the gravity waves at the center, but ar- 
gued that the waves would "break," i.e., become highly 
nonlinear, based on an estimate of the gravity wave am- 
plitude required for the waves to invert the stratification 
of the star. This argument is suppo r ted by the numer- 
ical si mulations of I Barker fc Ogilvid (|201Cll ) and I Barked 
(|20Tll) . 

Drawing o n Go odman fc Dickson's argument, 
lOgilvie fc LinI (|2007[ ) suggested that the dissipation of 
gravity waves may be much less efficient for planetary 
companions because the amplitude of the gravity waves 
are much smaller when they are excited by a planet 
rather than a star (given the factor of ^ 1000 mass 
difference between the two); in particular, they argued 
that the gravity waves would not overturn the stratifi- 
cat ion in the core and the effi cient dissipation invoked 
bv IGoodman fc Dicksonl (|1998[ ) would not apply in the 
case of hot Jupiter systems (see also I Barker fc Ogilvid 
120101) . In this paper, we show that nonlinear damping 
of the dynamical tide can be important even when 
the stratification of the star is not overturned. Linear 
theory thus does not apply even for solar-type stars with 
planetary companions. 

Searching for a more e fficient source of dampin g for 
tides in solar- t ype s tars, iWitte fc Savoniid (|2002j ) and 
lOgilvie fc LinI ()2007f) realized that far shorter length- 
scales, and hence enhanced damping, could be achieved 
by including the Coriolis force. Ogilvie fc Lin showed 
that Q ^ 10* can be achieved for solar-type stars based 
on the excitation and damping of short wavelength iner- 
tial waves in the convective envelope. This is, however, 
still not sufficient to explain the observed circularization 
of solar binaries. 

1.2. Structure of this Paper 

The remainder of this paper is organized as follows. In 
§ [5] we use perturbation theory to derive the nonlinear 
equations of motion and the key dimensionless coupling 
coefficients required to quantify the nonlinear couplings 
among stellar modes. Many of the details of these cal- 
culations are given in Appendix |Al The results in § [2] 
include two different ways of approaching the nonlinear 
problem. In the first f§ l2.1.ip . we expand quantities rel- 
ative to the unperturbed background state of the star. 
In the second (§ 12.1.2! and Appendix IA.6p . we expand 
quantities relative to the linear tidal solution. These two 
approaches are formally equivalent, but we will show that 
having both at our disposal is useful for understanding 
certain aspects of nonlinear tides in stars. In Appendix 
|B]we specialize to solar type stars and describe some of 
the key properties of the nonlinear coupling coefficients 



derived in § [5] and Appendix [Xj 

Having derived the underlying equations of motion, we 
then review the results of linear theory (§|3l) and summa- 
rize the key physical processes present in nonlinear tidal 
theory (§111). In § [5] we present the nonlinear stability 
analysis that is at the heart of this paper: this allows us 
to determine the conditions under which the linear tidal 
flow is nonlinearly unstable to the parametric instability. 
We also find a 'collective' parametric instability involv- 
ing many modes that grows significantly faster than the 
standard three-mode parametric instability. Section [6] 
presents examples of the types of nonlinear instabilities 
that afflict stellar tides, including both resonant nonlin- 
ear coupling (parametric instability) and nonlinear in- 
homogeneous driving. In §§ [7] and [8] we apply the re- 
sults of § [5] to solar type stars and assess the parametric 
stability of the dynamical and equilibrium tides, respec- 
tively. In § [9] we study off-resonant non-linear excitation 
of both stellar g-modes and p-modes by the equilibrium 
and dynamical tides. In §[10] we assess when the nonlin- 
ear interactions are sufficiently strong to invalidate the 
assumption that tides excite global normal modes. Fi- 
nally, in § [TT] we summarize our results, discuss their 
implications for solar binaries and solar-type stars with 
hot Jupiter companions, and discuss key directions for 
future research. 

The rigorous study of the synchronization and circu- 
larization of binary systems requires that one include the 
rotation of the star being tidally forced. However, includ- 
ing even the simplest case of rigid rotation greatly com- 
plicates the analysis. While the centrifugal force may be 
ignored to a good approximation, inclusion of the Cori- 
olis force changes the properties of waves with mode fre- 
quency smaller than rotation frequency, and new wave 
families appear. Specifically, eigenfunctions must be ex- 
panded in a sum over spherical harmonics, instead of 
a single harmonic for nonrotating stars. This greatly 
complicates the calculations and for simplicity we neglect 
Coriolis forces in this paper. 

2. STATEMENT OF THE PROBLEM 

We consider a primary star of mass M and radius 
R subject to a tidal acceleration from the secondary 
of mass M'. In a spherical coordinate system (r, 9, (f) 
centered on the primary, we take the orbit of the sec- 
ondary to be (Z?(i),7r/2, $(i)), where D{t) is the sepa- 
ration (neglecting backreaction on the orbit) and $(t) 
is the true anomaly corresponding to a Keplerian orbit 
with semi-major axis a, eccentricity e, and orbital period 
P^.b = 2TT[a^/G[M + M')]^/^ Defining the dynamical 
time of the primary Pq = 27r/a;o = 2tt{R^ /GM)^/^, the 
dimensionless strength of the tidal acceleration relative 
to internal gravity can be parameterized by 

M' fR 



M 



IQ- 



M' 



Pn 



/ ^orb 



(1) 



M + M'J\2.8hTj VlOdays, 

We treat the tidal acceleration VC/ ex £{GM / R^) as a 
small quantity compared to the internal gravity GM/R^ 
due to either the small mass or long orbital period of 
the secondary. Linear theory computes the response of 
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the star to the external tidal forcing to 0{e). In this 
work, we include the lowest order nonlinear terms, and 
hence consider effects at both 0{e) and 0{e^). As we will 
describe, the latter includes nonlinear driving of high- 
order modes by the dynamical tide and the equilibrium 
tide. 

2.1. Second order equations of motion 

To derive the second order equations of motion, let x' 
be the position of a fluid element in the perturbed star, 
X the position of the same fluid element in the back- 
ground state, and ^ the Lagrangian displacement vector. 
They are related by a;' = a; + ^{x,t). Likewise, p' is 
the true density and p is the background density. We 
write the internal forces due to pressure, buoyancy, and 
perturbed gravity that are linear in ^ as /i[^] and those 
due to leading-order nonlinear interactions as /2[^,^]- 
Explicit expressions for fi and /2 are given in Schenk 
et al. (2002; see their eq. [137]). 

The external forcing terms due to the companion can 
be deri ved from the inte raction term of the fluid Hamil- 
tonian ()Newcomblll962[ ). To second order, the fluid in- 
teraction Hamiltonian is 

Hint^ I d^x'p'{x')U{x') = I d^xp{x)U{x + ^) 



d^xp{x) 



u{x) + i-vu + -cev,v,u 



,(2) 



where d^x' p'{x') = d^xp{x) since mass is conserved. The 
first term in equation ^ can be ignored since it con- 
tains no dependence on ^. Taking a functional derivative 
SHint/S^ of this expression with respect to ^ leads to the 
tidal acceleration 



Itidc 



iSHint 
P 



S^ 



-VU-{^-V)VU. 



(3) 



The first term in equation ([3]) is the standard linear, time- 
dependent, inhomogeneous tidal forcing which acts to ex- 
cite oscillation modes. The second term is the nonlinear 
tidal forcing. Since it is linear in mode amplitude, it may 
lead to an exponential growth or damping of waves, as 
in th e Mathieu equation (see also iPapaloizou fc Pringld 

[mi). 

Gathering terms, the second-order equation of motion 
including linear forces, tidal forcing, and 3-wave nonlin- 
ear interactions, is then 



P^ = /iK] + /2[l,€]+patido 



(4) 



We consider two approaches to solving equation ^. 
Both approaches involve expanding the spatial depen- 
dence of all quantities in terms of the linear adiabatic 
eigenmodes of the star. In the first method, we expand 
quantities relative to the star's unperturbed background 
state. In the second method, we expand quantities rel- 
ative to the star's linearly perturbed state. As we will 
describe, each approach has its own conceptual and prac- 
tical advantages and we will make use of both throughout 
the paper. 

2.1.1. Method 1 



Following iSchenk et all (|200^ . 
dimensional phase space vector as 



expand the six- 



'^{x,t) 

i{x,t) 






qa{t) 



-lUJa$,a{x) 



(5) 



The eigenmode labeled "a" is specified by its frequency 
Wa, eigenfunction ^^(a;), and total amplitude qa- The 
sum over a runs over all mode quantum numbers, mode 
families, and frequency signs to allow both a mode and its 
complex conjugate. Plugging this expansion into equa- 
tion ^ and using the orthogonality of eigenmodes leads 
to a set of coupled oscillator equations for each mode (see 
Appendix [^ 



qa + iUJaqa = -JaQa + iuJaUa{t) 
6 



+lUJa} U:,,{t)qt + 



iuJaJ2'^*abcqWc-iQ) 
be 



The left hand side of equation ([6]) describes an uncoupled 
oscillator. The terms on the right hand side represent 
linear damping (7a), the linear (Ua) and nonlinear (Uab) 
tidal force, and three- wave coupling (Kate)- The latter 
two represent effects that come in at 0{e^). We choose 
a normalization in which, at unit amplitude, each mode 
has energy (see eq. |A2| ) 



Eo = 2u;i / d'xpC ■ ^a = 



GAP 
R 



(7) 



In terms of this normalizatioii. the coefficients in equa- 
tion dH) become (jSchenk et al.l(2002[ ) 



Uab{t) = -^J d^xpia ■ [i, ■ V) vf/, 



^abc 



(8) 

(9) 

(10) 



The tidal potential U is expressed in terms of spherical 
harmonics as 



U{x,t)-- 



E 

f>2.rj 



GM'Wirnr' 

D^+^{t) 



Yi, 



-im^(t) 



,(11) 



The Wim are 
— 2 harmonic, 



where Wim = 47r(2^ + l)-iy,„,(7r/2,0). 
nonzero only if £ — to is even; for the i 
which dominates for small R/D, W20 

W2±2 — \/3tt/10. Corrections to the tidal potential due 
to the extended mass distribution of the secondary oc- 
cur at quadrupole-quadrupole order (©(e"*) for a twin 
binary) and can therefore be neglected at C(e^). 

Plugging equation (fTTj) into equations dH) and ([9]) leads 
to the following dimensionless overlap integrals: 



and 



Jabtn. = j^i J d^xpia ■ iib ■ V) V {r'Y,, 



(12) 



(13) 
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In terms of these overlap integrals, the time-dependent 
coefficients in equation ([6]) can be written solely as func- 
tions of the orbit 






e.m 



£m 



Uab{t) = ^^ ^ Wl^Jab 



R 

W) 

R 



i*(t) 



(14) 



M 



Im 



D{t) 



i+i 



-im^{t) 



(15) 



For analytic work we find it convenient to expand the 
time-dependence of the Keplerian orbit as a sum of har- 
monic terms using the Hansen coefficients 



d) 



i+i 



— zm$ 



^ J2 xtiey 



ikQt 



(16) 



A:— — oo 



where 



Xr(e): 



2^ 



27r/n 

e 

"2 



dte 



ikQt — ini^ 



{3) 



e+1 



{£+1- 2m)4,m-i 



+{£+1 + 2to)4, 



m+1 



0{e^ 



(17) 



The linear tidal potential can be written in terms of the 



(fc)„-ifcOt 



Hansen coefficients as Ua{t) —J2k^<^ 



, where 



ini 



(18) 



and similarly for the nonlinear tidal coefficient Uab- 

The coupled mode amplitude equations can be inte- 
grated once the values of Wq, 7a, laem, JaUm and Kabc are 
known. In Appendix |A] we describe our calculation of 
each of these coefhcients. 

2.1.2. Method 2 

In method 1 we expressed the amplitude equation (eq. 
[5]) satisfied by each mode a in terms of the mode's total 
amplitude g^, i.e., the amplitude relative to the back- 
ground state of an unperturbed star. In method 2 we 
instead express the amplitude equation in terms of the 
nonlinear amplitude Va = qa — ga,iin, where qa.\m is the 
linear amplitude found by ignoring nonlinear coupling 
and the nonlinear tide (i.e., setting Kabc = Uab = in 
eq. [6]; see § |3] and eq. [25]). Although the form of 
the amplitude equation for r^ is somewhat less compact 
than that for g^, we find that this approach offers some 
advantages, particularly when studying the stability of 
the linear solution (see ij 15. 2p . 

Starting from the second-order equation of motion (eq. 
[4]), write the total displacement as a sum of the linear 
and nonlinear displacement ^ = ^y^^ -t- ^j^j. The linear 
displacement ^jjj^ is found by solving the hnear inhomo- 
geneous equations, as described in § lA.ll Noting that 
P^iin = /iKiin] - P^U, wc then have 



P^nl = /l Knl] + /2 [€li„: €li„] - 
+ /2Knb^„l] -p[(Clin 



-U)-V]VU. (19) 



Expanding ^j^j using equation ([S|) with (^,^, ga) -^ 
(^nb^nii^a) l^ads to au cquation for each mode's non- 
linear amplitude 

i-a + {iuJa + la)ra = i^a (V* + K^) 



:ib) ' b 



b 

+iWa^<,X^c*> (20) 



where 

Va{t)^~^J dPxpi,-{i,,,,-V)VU, (21) 

Ka{t)^^J d'xia ■ f, Kiin,^lin] , (22) 

Kab{t)^^Jd'xia-f,[^Un,^b]: (23) 

and Uab{t) and Kabc are given by equations ^ and (ITU)) . 
We describe our calculation of the coefficients Va, Ka and 
Kab in Appendix IA.6I 

2.2. Additional points 

Given the amplitude Qa of each mode, we can calculate 
the orbital evolution of the system (e.g., d and e). In Ap- 
pendix[C]we derive the expressions describing the orbital 
evolution assuming the excited modes are all standing 
waves. 

So far we have assumed that the tidally forced body 
is non-rotating. This may be a good approximation for 
tidal forcing of slowly rotating solar-type stars by close-in 
planets, however it is inadequate when studying circular- 
ization, where the body is pseudo-synchronized with the 
orbit. Within the approximations in this paper, where 
the Coriolis force is ignored, we can still use the cor- 
rect rotating frame forcing frequencies by replacing the 
stellar azimuth $ in the inertial frame with the rotating 
frame value ^rot = <& — mfls, where Qs is the rotation 
rate. This would alter the phase of the tidal potential 
to be exp[—i{kQ — mfls)t], where kQ — mfls is the co- 
rotating frame frequency. Note too that equation ([7]) is 
the mode energy for a non-rotating star, and it assumes 
that the potential and kinetic energies are equal. In a 
rotating star, equation ^ can still be used for normal- 
ization purposes, but to calculate the true energy the full 
kinetic and potential terms must be summed. Further- 
more, one must choose a frame (co-rotating or inertial) 
to evaluate the wave amplitude (g), orbital phase ($), 
eigenfunctions (£), frequency (w), and energy. 

Our treatment to this point has been for orbits with 
arbitrary eccentricity. Since the binary systems that are 
the focus of this paper typically have e ^ 1 and R <S^ a, 
for analytic work we use the dominant i — 2 component 
of the potential and harmonics up to linear order in e. 
For circularization the four harmonics are {m — 0,k — 
±1), (tti — ±2, fc = to), (to — ±2, fc — m — 1) and 
(to ~ ±2, fc = TO, + 1). For synchronous rotation fig = 57, 
the forcing frequency in the rotating frame is (fc — m)D, 
and the k = m term can be ignored since it has zero 
frequency. 

3. LINEAR TIDE 



WEINBERG, ARRAS, QUATAERT, & BURKART 



■10 



w 



w 



o 



■15 



-20 





A f\l = 9 




486 


- 


- 




\48 
\97 


J 




- 


- 






, 9'^A 


- 



0.01 



0.1 



1 



P« (days) 



10 



Fig. 1. — Linear energy E^ jin of individual modes (in units of 
Eq ; eq. [?] ) as a function of mode period Pa for the £ = k = m = 2 
harmonic of the tide in a non-rotating solar-type star with M' = 
M = Mq and Poih = 10 days. A line connecting the individual 
modes is drawn for clarity. The radial order n of the g-modes with 
period P ~ 0.1,0.5, 1,5, and 10 days are labeled. The broad peak 
centered on P ~ 0.05 days is the equilibrium tide and the narrow 
peak at P ~ 5 days is the dynamical tide. 



The linear response of the star can be found by ignoring 
the nonlinear coupHng (Kabc — 0) and the nonlinear tide 
[Uab = 0) in equation Q. Each mode then rings as an 
independent oscillator according to the equation 



whose steady-state solution is 



(24) 



/^^ Ua-kn-l-ya 

K— — C« 

The linear response can be broken up into a zero fre- 
quency equilibrium tide and a dynamical tide, qa,iin(i) = 

9a,oq(0 + '?a,dyn(i), where 



qa,.S)^ J2 Ui'h 



-~ikQ.t 



(26) 



and 



9a,dyn(i) = qa.\in{i) " qa,cq{t) 



= E 



UJa~kQ.- ija 



(27) 
^{k)^-rknt^ (28) 



The linear energy Ea^m — |<Za,iinP in an fa = 2 mode for 
the £ = k = m = 2 harmonic of the tide in a non-rotating 



solar-type star is 



Ea,lin — Ua 



ao" 



M' 



Pc 



orb 



Pa 

M + M'J V 10 days/ Vlday 



Eo 



-11/3 



(29) 



where A^ = Wq — kVl is the detuning and the second 
expression is valid for Pa ^ ^ day (i.e., when the linear 

overlap varies as laUm cc Pa ; see eq. |A16| ). In 
Figure [1] we show Ea^un for a solar binary with Porb = 
10 days. 

As we describe in § IA.21 the linear overlap integral 
laim, and thus Ua, is largest for low order £ = 2 modes 
with frequency uja ~ ^o- This is because the structure of 
these modes most closely resembles the shape of the tidal 
potential. For the orbits of interest, uto ^ H. and these 
low order modes are off-resonant with the tide; they are 
forced to oscillate at a frequency much lower than their 
natural frequency. These modes, seen as the broad hump 
in Figure [T] given by Ea.nn — Ea,eq = U"^, comprise the 
equilibrium tide and represent the nearly hydrostatic re- 
sponse of the star. Their net dissipation (by radiative 
damping in the stellar interior and/or turbulent damp- 
ing in the convective zone) yields solar b inary circulariza- 
tion times much longer than observed (jGoodman fc OhI 
[1991 . 

The dynamical tide is comprised of resonant high-order 
g-modes with uia — kfl (the peak at Porb/2 in Figure 
[1]). The dispersion relation of high order g-modes is 
(jja — ao{la/na) whcrc Qfo — 4 X lO"'^ rad s~^ for a so- 
lar model. The frequency spacing at fixed £a is therefore 
SuJa ~ u!'^/{ao£a) and statistically the detuning of the 
most resonant mode with the tide is \Aa/u}a\ ~ 10^"^ Pj^^ 
for £a = 2. We plot the linear energy of the most reso- 
nant mode for the £ — k = m = 2 harmonic of the tide 
in Figure [SI Numerically we find an average energy in 
the dynamical tide of 



E, 



dyn 



2 X 10^ 



-13 



M' 



M + M' 



Porh \ 

10 days/ 



-17/3 



Eo 
(30) 



(roughly the lower envelope of Ea,iin in Fig. [5]). For the 
periods of interest, this is ~ 100 times smaller than the 
linear energy of the low order modes with frequency ~ 
loq. However, despite their smaller energy, the dynamical 
tide in solar binaries is more dissipative than the equilib- 
rium tide for two reasons: (i) the resonant modes have a 
higher radiative damping rate owing to their short wave- 
lengths and, more importantly, (ii) the resonant modes 
appear to undergo such strong nonlinear damping in the 
core that within one travel ti me across the cor e they de- 
posit nearly all their energy (jGoodman fc Dick son 199^ 
see also § [TU]) . Nonetheless, the combined dissipation 
from the linear dynamical and equilibrium tides is ~ 3 
orders of magnitude too small to explain the observed 
circularization times of solar binaries. 
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4. NONLINEAR TIDE 

In linear tidal theory, the dynamics are dominated by 
^ = 2 modes since, by angular momentum conservation, 
only they can couple to the dominant ^ = 2 harmonic 
of the tide. In nonlinear theory, by contrast, modes 
with 1^2 can couple to the 1 = 2 harmonic of the 
tide. Angular momentum conservation requires only that 
the coupling coefhcient n^hc satisfy the selection rules 
1 4 - ic\ < ia < h + ic with £a + ib + ic evcu and 
TUa + rrib + rric = 0. The nonlinear overlap integral Jbcim 
is subject to the same selection rules but with £a and 
TTia replaced by the i and m of the tidal potential. This 
freedom in the i of the modes opens up a large region 
of the (£, n) parameter space that is inaccessible in lin- 
ear theory. In particular, since the g-mode dispersion 
relation is w ex i/n at the frequencies of interest, high 
i and n modes can be resonant with the tide and thus 
dynamically important. 

Nonlinear theory not only opens up the parameter 
space of modes accessible to the tide, it also modifies 
how modes interact with the orbit. This is because in 
nonlinear theory, the orbital evolution depends not only 
on the linear overlap laim but also on the nonlinear over- 
lap Jabim (see Appendix [C]) . For two modes a and b of 
similar (short) wavelength, Jabim 3> laim so that non- 
linear effects can be particularly important. Physically, 
this is because the coupled modes have an effective wave- 
length that is long and thus they induce a large density 
perturbation per unit mode energy. In this paper we fo- 
cus on determining the stability of the linear solution and 
do not attempt to calculate the rate of orbital evolution. 
Determining the full implications of the above effect is 
therefore deferred to a future paper. 

There are three types of nonlinear terms at C(£^), each 
the result of different forms of three- wave coupling. This 
is most easily seen in the amplitude equation for the non- 
linear amplitude r^ (eq. |20|): the three terms on the 
right hand side each represent a different form of three- 
wave couplingo In order of their appearance in equa- 
tion (1^ . the three types of three- wave coupling are: (1) 
"linear-linear coupling" (LLC) — the coupling of a nonlin- 
ear wave (or more precisely, the nonlinear correction to a 
wave) to a pair of linear waves (equilibrium tide plus dy- 
namical tide), (2) "linear- nonlinear coupling" (LNC) — 
the coupling of a nonlinear wave to a linear wave and an- 
other nonlinear wave, and (3) "nonlinear-nonlinear cou- 
pling" (NNC) — the coupling of a nonlinear wave to a 
pair of nonlinear waves. Conceptually, the terminology 
LLC and LNC is most useful when the linearly excited 
modes are near their linear energy. If nonlinear interac- 
tions drive these modes far from their linear energy, it no 
longer makes sense to call them linear waves and there 
is no longer a clear dichotomy between LLC, LNC, and 
NNC. We discuss this issue further in Jj^i 16.41 and [^ 

Like the linear inhomogeneous driving term Ua, the 
LLC term Va + Ka in equation ([^0]) acts as a (non- 
linear) inhomogeneous driving term. Thus, just as all 
modes that satisfy momentum conservation in linear the- 
ory have some linear amplitude ga.iin due to the driving 
term f/a, all modes that satisfy momentum conservation 

** These three terms can also be seen in the qa amplitude equation 
(eq. [6]) by substituting gj, = gj,^iin + rj, and qc = gc,lin + '"c into 
the nonlinear terms Y^bi'^ab + J2c i^abcQclqb- 



in nonlinear theory have some nonlinear amplitude ra 
due to the driving term Va + Ka- Insofar as the £ = 2 tide 
dominates, momentum conservation implies that only 
modes with £a = 0, 2, or 4 are driven significantly by 
LLC. 

The LNC term {Uab + '^Kab)rb has the form of a net- 
work of coupled Mathieu equations; it can thus lead to 
parametric driving. As with all parametrically driven 
systems, only modes that are sufficiently resonant with 
the driving frequency and have sufficiently small linear 
damping are unstable to LNC (see eq. [35]). This is in 
contrast to LLC, where all modes with nonzero couplings 
are driven regardless of their frequencies and damping 
rates. 

Finally, the NNC term Kabc^bfc has a quadratic depen- 
dence on the nonlinear mode amplitudes. Whereas the 
stability of LNC is independent of the amplitude of the 
nonlinear modes (it depends only on their frequency and 
damping rate) , only modes with sufficiently large nonlin- 
ear amplitudes are affected by NNC. NNC describes how 
the energy in nonlinear waves gets redistributed amongst 
other nonlinear waves and is thus related to the satura- 
tion of nonlinear instabilities. 

In this paper we focus on calculating the stability of 
the linear state. When the star is exactly in the linear 
state (i.e., the nonlinear amplitudes r^ are vanishingly 
small), the only nonlinear interaction is the LLC term 
Va + Ka- With just an infinitesimal perturbation from 
the linear state, the LNC term can become important, 
which represents parametric driving of 'daughter' waves 
by the 'parent' linear tide. In the next section, we carry 
out a stability analysis of LNC. Since our present focus 
is the stability of the linear state and not the saturation 
of possible nonlinear instabilities, we do not study NNC 
in this paper. 

5. PARAMETRIC STABILITY ANALYSIS 

In this section we derive the stability condition for 
parametric driving by the linear tide (the LNC term 
[Uab + 2Kab]rb of eq. [20]). The stability can be deter- 
mined in the usual manner: perturb the linear solution 
and see if the perturbations grow exponentially due to 
nonlinear forces faster than they damp. If they do then 
the linear state is unstable to parametric driving. In Ji l5.II 
we derive the criteria that determines the stability of the 
linear solution and in § 15.21 we determine the paramet- 
ric driving rate The on which the stability depends. We 
show that The can be decomposed into two distinct types 
of driving: driving by the equilibrium tide and driving 
by the dynamical tide (^^ 15.31 and 15.41 respectively; the 
stability of the dynamical tide and equilibrium tide are 
discussed in more detail in §§[7] and [8]). We also describe 
a collective form of parametric instability; in § 17.21 we 
show that as a result of this instability, the dynamical 
tide may drive a very large number of modes to large 
amplitude in a very short time (^ Porb)- 

5.1. Instability Criteria 

Consider the parametric (LNC) driving of a set of 
daughter modes {5, c} with ujh.c ~ i^/2 = kQ/2, where 
Lu is the tidal driving frequency. For such daughters, 
parametric driving will typically dominate over inhomo- 
geneous (LLC) driving, especially if ^b_c 7^ 0, 2, or 4. To a 
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good approximation we can therefore neglect the daugh- 
ters inhomogeneous driving term Vf, + K}, and, by equa- 
tion (|20l) . the amphtude equation for each daughter is 

h + [luJb + lb)n = luJb Y, put) + 2KUt)] rl. (31) 

c 

Consider a particular harmonic uj of the tide such that 
Ubc{t) = Ubce~^'^'^ and similarly for Kbc{t). If we plug 
in rb{t) = (5b(i)e*"*/^ and similarly for r^ then the har- 
monic time dependences cancel, yielding 



{iAb + 7b)Qb^»y^ 



1/2 



HcQ:, (32) 



where we defined the daughter pair "driving rate" 

The = y/i^b^c [Ubc + SiVTfcc) . (33) 

Writing these equations as Q = -ffQ, the solutions are 
Q oc exp(si) where the eigenvalues s are the solutions 
of the characteristic equation det(_ff — si) — 0. The 
system is unstable if there is an s such that Re{s) > 0. 
For N coupled daughter modes, _ff is a non-symmetric 
2N X 2N matrix with 2x2 block components Hjk — 

{ijjj/uJk)^''^Ajk + 6jkBk, where the indices j,k run over 
the daughter modes, and 



ijk 



■Im(r,fc) Re{Tjk) 

Rc(r,fe) -im(r,fc) 



Bk 



-Ik Afc 
-Afc -7fc 



(34) 



Here A^ — Uk+uj /2 is the daughter detuning rather than 
the previously defined parent detuning Aa = Wq — w. 

While in general the eigenvalues of the characteristic 
equation must be solved for numerically, we give exam- 
ples below where the eigenvalue expressions are simple 
and offer insight into how the stability depends on iV, 
Wb, 7b, Ab, and Fbc- From these examples and numerical 
experiments we deduce an approximate stability criteria: 
a set of A'' daughters with Wbi^c > is unstable if each 
pair (&, c) in the set satisfies 



N\Tbc\ > / 



7ft7c 






(35) 



where Abe = Ab + Ac = u + ujb + ^c- When sufficiently 
far from the stability boundary, the growth rate of an 
unstable collection of modes is ~ A^|Fbc|- If WbWc < 0, the 
system can only be unstable if Fbc is asymmetric inb ^ c. 
However, since Ubc and Kbc are symmetric in 6 o c, so 
is Fbc- High-frequency modes with ujbUJc < and small 
Abe (i-e., a pair with a beat frequency \ujb + ud — \i^\) 
are therefore stable. 

We emphasize that to be collectively unstable, and thus 
have a growth rate ~ A^jFbd, each pair in the set of N 
daughters must satisfy the inequality (1351) . The inequal- 
ity's dependence on N shows that the daughters can be 
collectively unstable even if each daughter pair would be 
stable on its own (i.e., even if |Fbc| < \/A^^ -|- 7b 7c for ev- 
ery pair) . Conversely, a subset of pairs in the set can be 
unstable even if the set as a whole is collectively stable. 
This implies that even if the tide excites A^ ^ 1 daugh- 
ters, these daughters are not necessarily undergoing a 
collectively instability. Rather, they may each be under- 
going the standard parametric instability with a growth 



rate <^ A^|Fbc| (if the daughters are coupled to each other 
in a chain, they may all grow at a single, coherent rate 
<^ A'^IFbcl corresponding to one of the eigenvalues Re{s)). 

The case A^ = 2 is the one most often found in the 
literature ("three- wave coupling"). As we show in §§ [6] 
and [71 the collectively unstable case with N ^ 2 results 
in mode dynamics that are very different from the N — 2 
case; in particular, the driving rate of the N daughters 
is much more rapid (by a factor of ~ N/2). In § 15.41 we 
show that due to collective driving, the global growth 
rate of daughter modes approaches their maximum local 
growth (in the case of dynamical tide driving, the latter 
is the local growth rate within the small driving region 
near the center of the star). 

Physically, collective driving occurs because the A^ un- 
stable daughters each have comparable group velocities 
and thus their superposition within the growth region re- 
mains coherent over a growth time. They therefore act 
as a single "mode" that is much more strongly peaked 
than the individual eigenmodes. This results in an ef- 
fective coupling coefficient ^ A^ times larger than the 
three- mode coupling coefficient Kabc- 

The astrophysics literature has traditionally focused 
on three - wave coupling u sing global normal modes (how- 
ever see lRvu et al.1119961 ). What we call collective driv- 
ing is not necessarily captured by such an analysis, even 
though it is simply a consequence of spatially localized 
coupling of many modes. In particular, the rapid 'col- 
lective' growth at nearly the maximum parent shear rate 
(see § 15.41 below) is not captured in standard calcula- 
tions of three mode coupling coefficients and growth rates 
(e.g.. IKumar fc Goodmal]ll99l Iwlfc GoldreichI l200l 
[Weinb erg & Q uataertI 120081 )^ The reason is fundamen- 
tally that the standard basis set of global stellar nor- 
mal modes does not necessarily capture the most rapidly 
growing daughters in the system, which are in fact a 
superposition of standard global stellar normal modes 
(those that are strongly peaked near where the parent's 
shear peaks). In future work, it will be important to 
revisit previous astrophysical applications of parametric 
resonance given this limitation of previous work. 

The stability analysis of equation (|31|) applies only if 
the nonlinear interactions are sufficiently weak that the 
tide excites standing waves (i.e., global normal modes). 
As we describe in § llOi if the nonlinear interactions are 
so strong that the daughters' local growth rate (eq. [62] ) 
is larger than the inverse of their local group travel time 
(eq. [61]), then the driving must be treated as a local 
interaction in which the tide excites traveling waves. We 
note, however, that the instability criteria (eq. [35]) is 
independent of the uncertainty of traveling versus stand- 
ing waves. This is because the daughter growth rates are 
infinitesimal at threshold; the growth times are therefore 
always much longer than the group travel time very near 
threshold (a system is near threshold for certain values of 
the companion mass and orbital separation; see Fig. [7]). 
Thus, while collective driving may lead to growth rates 
that are much faster than that of standard three-mode 
coupling, a collectively unstable system is not necessar- 
ily in the regime in which the excitation of the daughters 
must be treated using traveling waves rather than global 
normal modes. 

5.1.1. Examples 
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(i) N daughters, 7b = 7, Afc = A, The = F : Two of the 

eigenvalues are s = —7 ± -^/lA^rp — A^ and the 
other 2{N — 1) are s = —7 + iA. The system is 

therefore unstable if A^|r| > ^7^ + A^. We obtain 
the same result if instead Tjk ~ (— l)-'+'^r; the sign 
of Kajk varies in such a manner for daughters that 
are An = 1 neighbors. 

(u) N = 2 daughters, Tbc = T, Tbb ^ Tec = : This 

system is often considered in the literature (e.g . , 
IWu fc Goldreichl[200l IKumar fc GoodmanllTOOeh . 
The four eigenvalues are 



s = -<i -(76 + 7c)±i(Ab- Ac)± 



4|rp-AL 



+ (76 - Icf ± 2iAhc(7c - 7h) 



1/2 



(36) 



with the first and third "±" linked. The system is 
unstable if 



m>V^mc 



1 



A? 

DC 



{lb + Icf 



1/2 



(37) 



Since this case does not include self-coupling, we 
can have lS.bc <C Ab, Ac. Thus, even if the daugh- 
ters' frequencies are not individually close to — w/2, 
their sum can be close to —uj and the pair can have 
a low parametric threshold. 

5.2. Parametric Driving Rate 

We now discuss the parametric driving rate of a daugh- 
ter pair (6, c). In this case N = 2 and the driving rate 
A^lLfccI ~ iLfjcl (eq. [33]). Expanding the coefficients 
in terms of the tidal harmonics, the driving by the fc'th 
harmonic is given by 



.(fc) 

- be 



U, 



(fc) 



- y/oJb^. 
Im 



2K 



(fc) 



R 



-(eq) 



-(dyn) 



,(38) 



where we define the equilibrium tide and dynamical tide 
driving rates 



=;(eq) 



..(<=q) 



-(dyn) 



•'be = i^bc V^b^c, •'be = "-be 

and their dimensionless coupling coefficients 



Ac^f^^V^J^- (39) 



,(cq) 



'Jbe 



2k 



(I) 



2k 



{H.cq) 



(40) 

(41) 

In Appendix [B] we describe the properties of the coeffi- 
cients Kli and Ki/ for the case of a solar-type star. 



(dyn) _r, (ff,dyn) 
'^be —'^'^be 



^bc 



^be 



Here Ki„ and kL ' ^'^' are the homogeneous parts of 



.(I) 



"be """^ "-be 

the linear tide coupling coefficient an d k^J is the inho- 
mogeneous part. As we explain in § lA.Sl the homoge- 
neous parts are found by directly replacing mode a in 
our final expression for Kabe (lines IA55IIA62)) with ^ 
and Cdyn, respectively (this expression for Kabc is simi- 
lar to that typically used for three-mode coupling in the 
literature, with a few corrections). The inhomogeneous 



part k\J corresponds to terms in the coupling coefficient 
that are not present in the existing treatments of three- 
mode coupling in the literature. They arise because we 
are considering non-linear coupling among tidally foreed 
modes, rather than freely oscillating modes. More specif- 
ically, the inhomogeneous terms arise because in deriving 
our final expression for Kaho we assume all three modes 
satisfy the homogeneous equations of motion (eqs. |A33[ 
IA34| '). Here, however, mode a is the linear tide ^ij„ which 
instea d satisfies t he inhomogeneous equations of motion 
(eqs. |A33| IA34| with S(j) ^ Sc/) + U); substituting the 
equations of motion which include the inhomogeneous 



be 



term U introduces additional terms which we call k 

Because kj^J does not contain an explicit dependence on 
^ijjj, it does not separate into an equilibrium tide piece 
and a dynamical tide piece. 



(cq)| 



.(dyn) 



We find that for a solar-type star |r^7^| < |r 
However, the driving rates are not the only diagnostic 
of an instability's importance. In particular, since the 
majority of the linear tidal (interaction) energy is stored 
in the equilibrium tide (§ [3|), an instability of the equi- 
librium tide could in principle lead to faster orbital evo- 
lution than an instability of the dynamical tide. Equi- 
librium tide driving thus offers a potentially important 
source of energy loss that is distinct from the dy i iamic al 
tide driving considered by iGoodman fc DicksonI ()1998[ ) . 

5.2.1. Relation between global and local driving rates 

While the exact expressions for the coupling coeffi- 
cients are complicated, intuitively one expects the local 
driving rate of short wavelength daughter waves to ap- 
proximately equal the driving frequency oj times the local 
'shear' of the parent wave ~ d^r/dr. Indeed, in ? IA.4.1I 
we show that the dominant terms in the full expressions 
for K,abb are approximately equal to 



dn. 



M 



dlnr 



T 



dEb d£,rM 
d In r dr 



(42) 



where T is an angular integral (|T| ~ 1) and dEb/dr ~ 
pr^N'^S.rb is the radial energy density of the daughter 
waves. Assuming high-order daughters (eqs. |A3II A5| ) . 
this implies that the global driving rate of a self-coupled 
daughter 5 by a parent a is approximately given by 



. bb = LOKabb 



-(local) 



iKlT, 



(local) 
bb 



dr 



JNT 



—(local) 



bb 



dlni 



I^r.bdr 



J Ndlni 



.(43) 



where F^^ (r) = ujTd^r,aldr is the local daughter driv- 
ing rate, Vr,b — ^b/kr,b is the daughter group velocity, 
and the limits of integration are the inner (ri ) and outer 
(r2) turning points of the daughter (note that TV here is 
the Brunt- Vaisala frequency, not the number of modes). 
The global driving rate is therefore the average of the lo- 
cal driving rate weighted by the time the daughters spend 
at each radius (or, equivalently, weighted by N/r). 

5.3. Equilibrium tide driving rate 

Naively, the nonlinear tide Ubc and the equilibrium 
tide part of the three wave coupling ^^ KabcQa.cq have 
a completely different origin. However, in Appendix lA. 51 
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we show that the leading order terms cancel in the in- 
tegrand of Jbcim + 2k[P (in eq. 



for nj^^ )■ Specifi- 
cally, the fractional difference between Jhcim and SkJ^^ is 

~ (lo/NY ~ 10~^^io^ ^^'^ there is thus a large cancella- 
tion between these terms in the growth rate of daughter 

pairs due to the equilibrium tide, Fj,^ . Physically, this 
is because nonlinear tidal driving {Uhc) and internal non- 
linear driving by the equilibrium tide (^^ Kabc9a, oq) are 
fundamentally part of the same process; together they 
describe the nonlinear driving of daughter modes by the 
nearly hydrostatic response of the star to its companion. 
In their paper on the parametric excitation of mo des in 
close binary systems, iPapaloizou fc Pringld ()1981l ) con- 
sidered driving by the nonlinear tide Ubc- However, they 
neglected three wave coupling ^^ Habcla- Their analysis 
therefore overestimates the equilibrium tide driving rate 
by a factor of ^ 10^ for solar-type stars. 

In Appendix IA.4.21 we show that dK\^^ /d In r is ap- 
proximately given by equation (|42p . where now mode 
a represents the equilibrium tide and d^r,a/dr « 
A(r/_R)^+^ at r ~ i?. The magnitude of d{Ji,ctm + 

2k\ J) I d\ \Yr is easily shown to be of the same order (see 
eq. |A71| '). Thus, the local daughter driving rate by the 

equilibrium tide is F^f, ^ uj{r / RY^^ . Although the 
local driving rate is ~ oj at r ~ i?, the global driving 

rate rf,^ (eq. [43]) is considerably smaller than this 
local value. This is because the global driving rate is 
weighted by the time the daughters spend at each radius 
and the daughters propagate more slowly at small radii 
(ur,b c>c r/N), where the shear is small (equivalently, the 
energy in the daughter modes peaks at small radii, where 
the shear is small). For a solar binary, the local driving 
rate by the equilibrium tide is everywhere much smaller 
than the inverse of the daughters' group travel time (see 
the next paragraph). It is therefore appropriate to treat 
the daughter excitation using global normal modes and 
their growth rate is given by the global driving rate. 

We can obtain an approximate expression for Ff,t, 
by noting that for a solar model N « 100a;o(r'/i?) for 
r < 0.05i? and N « Swq (to within a factor of two) be- 
tween r « 0.05i? and the convection zone re — 0.7R. 

With these approximations we find F;,^ « 0.2T{i + 
^)^^{rc/RY^^oj « O.Olcj, in good agreement with the 

full integration of k^j^ given in § IB.2R The global driv- 
ing rate by the equilibrium tide is therefore a factor of 
~ 100 slower than the local driving rate at r ~ i?. For 
the dominant harmonic of a synchronized binary (see eq. 



]), 



-^(oq) 
bb 



eel bb 



10" 



eM' 

M + M' 



Pfo'yr-S (44) 



where e = {M' /M){R/aY+^ (eq. [T]). 
We can use this estimate of F^^'^ to evaluate the sta- 

^ The integral in the numerator of equation 1143 ll is dominated by 
the contribution at r ^ ij, where to a good approximation g oc r~-^ 



and thus ^r,oq — 
g oc r and ^r.cq 
integral. 



-U/g (X r^+2 (gee eq. |A12j '). The core, where 
r^~^, contributes only a few percent to the 



bility of the equilibrium tide. The linear damping rate of 
a high order daughter is jb ~ 'T-b/^KH, where for a solar 
model ixH ~ 3x 10^ yr and rib — 500^f,Pio.b. If we ignore 
the daughter detuning we find that the equilibrium tide 

is unstable (F^^'' > 76) for orbital periods 



P < 4£, 



-2/5 



eM' 



M + M' 



1/5 



days. 



(45) 



While there may be brief intervals during which the de- 
tuning Abe ^ lb and the above limit applies, the detun- 
ing, on average, sets a more stringent constraint. This is 
because the equilibrium tide driving rate peaks strongly 



for self-coupled pairs (F^^ <C ^bb ^^^ \^b 



(eq) 



^(eq) 



see t; IB.2l) : thus, although there may exist daughter pairs 
with very small detuning Abe, in general such pairs have 
\nb — nc\ ^ 1 and they couple very weakly to the equi- 
librium tide. Using the relation for 7b and self-coupled 
detuning Abb ~ '^b/2?T.b ~ 0.2^^^ Pj^ ^ yi'~^ and minimiz- 
ing the right hand side of equation ((35)) with respect to 
Ub, we find that, on average, the equilibrium tide is only 
unstable to parametric resonance for 



P < 0.5 



eM' 
M + M' 



3/7 



days. 



(46) 



These estimates agree well with the more exact treatment 
given in § [S] 

Although we defined the parametric growth rate Fbc 
in terms of the method 2 coupling coefficient Kab 
(S I2.1.2p . we could have equivalently defined it in terms 
of the method 1 coefficient Kabc (§ I2.1.ip . To do so, 
write the coupling coefficient in functional form Kabc ~ 
filial ib, id- Then 



Kbc^K[^lin,^b,^c\ == '^ 



/ , "ZaJinSa' Sbi Sc 



= X! '* [^"i' ^fc' €c] QaMn = ^ KabcQaMn- (47) 

a a 

Our definition of Fbc (eq. [33] ) would thus involve a sum 
over all parents, including the low-order modes that com- 
prise the equilibrium tide and the high-order modes that 
comprise the dynamical tide (see §[3]). In principle, we 
could calculate Fbc by numerically summing the product 
of each parent's Kabc and qa.nn (eqs. [I0l[25]). However, 
in practice our calculations of individual Kabc are only ac- 
curate to a part in ^ 10"^ (even after c onsid erable effort 
to minimize the numerical error; see ij IA.4|) . Since Ubc 
and the equilibrium tide part of 2Kbc {— 2 J2a '*abc9a, cq; 
see eq. [26]) cancel to a part in ^ 10^, we cannot ac- 
curately calculate the equilibrium tide part of Fbc by a 
term-by-term sum over modesjj We find that the only 
way to accurately calculate this part of Fbc is to use our 
method 2 formulation, in which we use the full solution 
to the inhomogeneous linear equations ^jj^ rather than 
its mode decomposition ^^^^ = ^ 9a,iin^a- 

5.4. Dynamical tide driving rate 

^ We have verified that they cancel to the precision of our sum 
over modes. 
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Unlike the shear of the equihbrium tide, which peaks 
at r ~ i?, the shear of dynamical tide peaks in the core 
(since dS,r/dr ~ fc^^r (x r~^). The driving of daughters 
by the dynamical tide therefore occurs primarily near 
the dynamical tide's inner turning point ri. If we de- 
fine the nonlinearity parameter of the dynamical tide 
A = max(fcr^r) PC rr^ fsee iGoodman fc Dicksonlll998l : 



MU 



lOgilviefcLiiil 120071 ). then by equation (gll), the global 

driving rate of daughters is rj,^''" ~ wTAri/r2, where 
r2 ~ re — 0.7R. The inner turning point is located 
where oj ~ N ~ lOOojor/R and thus ri/r2 « 0.2/n « 
2 X 10~^P^Q^, where n is the radial order of the dynami- 
cal tide_(thisjm£liesnj^^ oc LuA/n, a scaling also noted 
by [Barker fc Ogilvii 120 lit ) . The global driving rate of 
the dynamical tide is therefore 



n(dyn) 
bb 



0.01 



fM'/M\ p_n/6 
^ 10-- 



10 



yr 



(48) 



jl/6 



where we took A ~ 470(M'/M)P]^q" corresponding to a 
tide raised i n a slowly r otating solar-type star by a planet 
of mass M' (jO gilvie fc Lin 2007). The ratio ri/r2 is ap- 
proximately the ratio of the global driving rate to the 
maximum local driving rate (for comparison, the global 
driving rate by the equilibrium tide is ^ 1% of the max- 
imum local driving rate). The above estimate of T^^''" 
agrees well with the more exact value found by multiply- 

ing our numerical integration of Kabc with -E^' (see also 

ea s. | A67I and [30]); it also agrees with the rate found 
bvlBarkcr fc Ogilvii (j2011h . At P = 10 days, the driving 
rate by the dynamical tide is ~ 10^ times faster than the 
driving rate by the equilibrium tide (eq. lU). 

The above estimate assumes that only a single daugh- 
ter pair is excited and thus ignores the possibility of col- 
lective driving. The collective driving rate due to the 

(dyn) 



dynamical tide is « A^r^^ . In Appendix IB. II we show 

that A^ « n (see also §[72]). Thus, ^F^^''"^ « 0.2ujTA, 
i.e., the collective global driving rate is approximately 
equal to the maximum local driving rate of individual 
daughter pairs r^^°'^^ ,max) _ uXj\^ 

6. ILLUSTRATION OF NONLINEAR 
INSTABILITIES 

In this section we consider simple coupled networks in 
order to illustrate how nonlinear interactions can redis- 
tribute the energy of excited modes. In the first three 
subsections we give examples of parametrically unstable 
systems (i.e., LNC; see §[4]) and in the last subsection we 
give an example of inhomogeneous driving (LLC). More 
specifically, in § 16.11 we consider a parametrically unsta- 
ble three-wave system consisting of a linearly resonant 
parent a coupled to a pair of daughter modes 6, c with 
frequency 2:^ uJa/2. We assume that this is the only source 
of daughter driving (i.e., we set Utc = 0). Because the 
dynamical tide is comprised of linearly resonant parents 
(§[3]), this example helps illustrate the parametric insta- 
bility of the dynamical tide discussed in § [7] In § 16.21 we 
consider the same three- wave system as in § 16.11 but in- 
clude nonlinear tidal driving of the daughters {Ubc y^ 0). 
Because the equilibrium tide drives daughters at a rate 
that depends on an effective Utc (§ 15. 2p . this example 
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Fig. 2. — Effect of three- wave coupling on a linearly resonant g- 
mode without nonlinear tidal driving {Uab = 0). The parent (solid 
black line) and daughter modes (dashed gray line; since they are 
similar, only one of the dau ghters is shown) are initially at the 
linear tide solution (eq. |29|). For the parameters of the system 

1/2 _o 

we use ia = £i, = ic = 2, V^c = f^abc^^^'lin^ = 10 ^fJ and ■ya = 

"fb = Ic = Ajjc = 10~*n. These parameters are similar to those of 
a solar binary. The linear solution is unstable and the daughters 
initially grow on a timescale of ~ l/F^j. ~ 100 orbits. After ~ 10* 
orbits the system reaches a new steady state nonlinear equilibrium 
(dotted lines; see Appendix [DJ. 



helps illustrate the parametric instability of the equilib- 
rium tide discussed in § [8] In § 16.31 we show an example 
of a collectively unstable system in which a single lin- 
early resonant parent is coupled to A'^ daughters all with 
frequency ~ Wq/2. We show that the collective instabil- 
ity can lead to the very rapid growth of a large number 
of daughters and thereby dominate the mode dynamics. 
Finally, in § 16.41 we consider a linearly resonant parent 
a coupled to itself and a pair of daughter waves with 
frequency ^ uJa- This example helps illustrate nonlin- 
ear inhomogeneous driving by the linear tide discussed 
in §§[land[i 

6.1. Three-wave parametric instability 

The three- wave coupling considered here involves a lin- 
early resonant parent mode a coupled to two linearly 
driven daughter modes 5, c of approximately half the par- 
ent's frequency (this is sometimes re ferred to as th e para- 
metric subharmonic instability; see iMiiller et al.l [l986) . 
We artificially set the nonlinear tide coefficient to zero 
for all three modes {Uij = for {i,j} = {a, 6, c}). The 
amplitude equation for the parent is thus 

Qa + {i^a + la)qa='i^a [Ua{t) + 2Kabcqlq*c\ (49) 



and similarly for the daughters. We see by equations (|33]) 
and ([37]l that the daughters are unstable if the parent's 
linear energy Ea^iin = |Qa,iinP is above the parametric 
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chose the ratio of their driving rates to be 1/4, setting 



logj^Q Time (orbits) 



Fig. 3. — Same three mode system as Figure[2]but now including 
self-coupled nonlinear tidal driving of the daughters. We use Utb = 

Ucc = '^abc-E^ iin/4 = 2.5 X 10"'^, an artificially large value for a 
solar binary. For comparison, the thin dotted lines shovf the parent 
and daughter of Figure [2] in which Utb = Ucc = Ui,c = 0. The 
system is unstable and parametric coupling to the parent cannot 
saturate the exponential growth. In reality, the equilibrium tide is 
only unstable to parametric resonance in solar binaries with Porb ^ 
2-5 days {§[8ll. 

threshold 



E, 



th^ 



7&7c 



4'«a6c^fc^c 



A? 
be 



:10 



-20 



10 days 



(7b + Icf 

2 
1 



A? 

be 



iib + icY 



Eo (50) 



where in the second hne we assumed Pb — Pc — Ppxh and 
plugged in the analytic expressions (IA14[) and (IA67P for 
7f,,c and Kabc- 

In Figure [H we show the evolution of such a para- 
metrically unstable three mode system. The modes are 
initialized at their linear energy (eq. [15] )• The lin- 
ear solution is unstable and the daughters grow at a 

rate Vbc ~ '^abcEjy^^^- As we show in Appendix [Dl 
this system eventually settles into a stable equilibrium 
with energies approximately given by Ea — £'thi Eb — 

{'^ciObhbi^cY^'^\Ual'2.K.abc\-, and Ec = {jbUJc/jcUJb)Eb. 

6.2. Nonlinear tidal driving 

Just as modes with frequencies ~ uj are linearly reso- 
nant with the tide through the Ua term in equation ([6]), 
modes with frequencies ~ uj/2 are nonlinearly resonant 
with the tide through the Uab term. The daughters in 
the previous subsection are therefore subject to nonlin- 
ear tidal driving in addition to three- wave coupling with 
the linearly resonant parents. 

In Figure [3] we show the evolution of the same three 
mode system as Figure [5] but now including nonlinear 
tidal driving of the daughters. Such driving is analogous 
to driving by the equilibrium tide if it were to be para- 
metrically unstable. In order to contrast the driving by 
the nonlinear tide with that by three wave coupling, we 



Ubc — Ubb 



Uc, 



^abcE]!lj^ = 2.5 X 10-3 in the 



daughters' amplitude equation. Our choice of Jbcim in 
Figure [3] is a factor of ~ 10^ larger than the true ef- 
fective value of the equilibrium tide coupling coefficient 

*^bc ^'-'^ solar binaries (see ij 15.21 and eq. [10]). We use 
an artificially large value in order to illustrate the nature 
of the nonlinear instability. 

1/2 

Since Ubc < i^abcE^Yuv tti6 daughters' driving is ini- 
tially dominated by three- wave coupling to the parent. 
There occurs a short-lived nonlinear equilibrium with en- 
ergies similar to the equilibrium of Figure dj However, 
this equilibrium is destabilized by Ubc and after w 300 or- 
bits the nonlinear driving by Ubc begins to take over and 
the daughters grow exponentially at a rate Vbc — Ubb^- 
Three- wave coupling to the parent does not saturate the 
daughters; instead the daughters drag the parent with 
them and the parent also grows exponentially. 

6.3. Collective instability 

In Figure [4] we show a system of A^ = 100 daughters 
that are each coupled to each other and the same linearly 
driven parent as in Figure [5] (a total of 10^ couplings). 
The magnitude of the daughter damping rates and de- 
tunings are 7, A^c — 0.1J7, corresponding to high i, some- 
what resonant daughters. We find that the stability of the 
daughters is not sensitive to their initial conditions. In 
order to illustrate this, we assign each daughter a ran- 
dom initial phase and amplitude such that their initial 

energies vary over 10 orders of magnitude. 
1/2 
Since Vbc = KabcE^y^^^Vl = O.Olfi, on their own each 

daughter pair would be stable according to the three- 
wave criteria (eq. [37]). However, since the set of 
N daughters satisfies the instability criteria of equa- 
tion (|55]) . they are collectively unstable. The daughters 
with the lowest initial energy grow fastest and within a 
few orbits the entire set is growing coherently at a rate 
A^lEfjcl ~ ^- This is a factor of ~ A^ times faster than 
that of the standard three-wave parametric system con- 
sidered in § 16.11 Once coherent, the daughters oscillate 
in unison at a frequency a;/2 and eventually settle into 
a stable equilibrium with energies that are all equal to 
within w 30%. 

The rapid daughter growth rate implies that the lin- 
early driven parent loses energy at a much faster rate 
than the parent in the three- wave system of § 16.11 This, 
in turn, implies that a collectively unstable system can 
potentially drain energy out of the orbit much faster. As 
we discuss in § 1101 this may imply that the driving is in 
the traveling wave (rather than standing wave) limit for 
much lower mass companions than indicated by ordinary 
three- wave coupling. In ij 17.21 we show that the dynami- 
cal tide due to Jupiter-mass planets in few day orbits is 
collectively unstable, driving A^ sa lO^Pio daughters with 
I K, 3OP20 to significant amplitudes within a few orbits. 

6.4. Three-wave inhomogeneous driving 

The examples in the last three subsections all illustrate 
the excitation of modes through the parametric instabil- 
ity (LNC). As described in § |H the linear tide can also 
excite modes through nonlinear inhomogeneous driving 
(LLC). In terms of the rb version of the amplitude equa- 
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Fig. 4. — Collective parametric driving of a set of A'^ = 100 
coupled daughters. The parent is the same as in Figure [2] with 
Tbc = O.OIO while each daughter's damping rate and detuning is 
increased by a factor of ~ lO'^ to ^i, — Ic — ^bc — O.lfi. The 
system is collectively unstable even though each daughter pa ir is 
stable according to the three-wave instabiUty criteria (eq. 1371 ). 



tion (eq. [20]), nonlinear inhomogeneous driving of a 
mode b arises from the coupling coefficient Vb + Ki,. This 
coefficient accounts for the full linear tide and thus inher- 
ently includes a sum over all linear parent modes. Fol- 
lowing the simple examples of the previous subsections, 
here we instead consider inhomogeneous driving by just 
a single linearly resonant parent mode a. We couple this 
parent to itself and two high-frequency, non-resonant, 
daughter modes 5, c. In terms of the qa version of the 
the amplitude equation (eq. [6]), we have 

Qa + {iuJa + la)qa = i^a [Ua{t) + 2 (Kaabqt + KaacQc) 9q] 
qb + {iUJb + lb)qb = i^b [Ub{t) + Kaab{q*a)'^] 
qc + {iuJc + 7c)<?c = JWc [Uc{t) + Kaaciql)"^] ■ (51) 

We show an example of such a system in Figure [S] Un- 
like the parametric instability considered in § 16. 1[ which 
is only unstable if Ea,iin > Eth (eq. [50]), the linear solu- 
tion of this system is always invalid because the parent 
appears as an inhomogeneous driving term in the daugh- 
ter amplitude equations. Thus, just as the steady state 
solution to the linear equation for high-frequency modes 
driven at a frequency lu <ti uja is qb^ua ~ UbC^'^* (eq. [55] 
assuming jb ^ ^b), the daughter steady state solution 
here is ^6 w KaabEac'^'''^* (assuming KaabEa > Ub). For 



■12 



Figure [5] we chose coupling coefficients 



e: 



-1/2. 



the parent's steady state energy is therefore only 
slightly smaller than -Eajin- Depending on the tidal 
factor £, the actual inhomogeneous driving coefficient 
Vb + Kb can be so large that the parent's steady state 
energy is very different from -Eajin (see § [9]). In that 
case, the simple view of "linear-like" driving no longer 
holds and the other forms of three-wave coupling (LNC 
and NNC) can become important. 
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Fig. 5. — Nonhnear inhomogeneous driving of two high-frequency 
daughter waves by the same Hnearly resonant parent as in Figure 
[2] The daughters have frequencies wj, = WuJa and uJc = lOOwa 
{dashed and dashed- dotted lines, respectively) and 7a = 7(, = 7c = 
10~*n. The daughters couple to the self-coupled parent (solid 
line) with cou plin g coefficients K^ab = i^aac = 10^. The linear tide 
solution (eq. 1291 : dotted line shows E^ un) is used as the initial 
condition. The system settles into a nonlinear equilibrium in which 
the daughter energies are much larger than their linear values while 
the parent's energy is only slightly smaller than its linear value. 



7. PARAMETRIC INSTABILITY OF THE 
DYNAMICAL TIDE 

As discussed in § [S] the linear tide can be decomposed 
into a dynamical tide (high-order modes resonant with 
the tide) and an equilibrium tide (low-order, off-resonant 
modes with large overlap integrals I aim)- In this section 
we show that for the close binary systems discussed in §[1] 
the dynamical tide is parametrically unstable to nonlin- 
ear three-wave interactions. We first (? l7.ip assume that 
each daughter is coupled to the dynamical tide and only 
one other daughter {N < 2) and determine the param- 
eter space over which the daughters are unstable. We 
then (§ 17. 2p allow each daughter to couple to the dynam- 
ical tide and many other daughters (TV ^ 2). We show 
that these collections of modes are collectively unstable 
and have growth rates that can be orders of magnitude 
larger than the TV < 2 systems. 

Since the dynamical tide is comprised of several par- 
ents on either side of the linear resonance, in principle, 
one should allow each daughter pair to couple to more 
than one parent. However, we find that the most linearly 
resonant parent dominates the dynamics and for simplic- 
ity we only show results for driving by this single parent. 
Furthermore, although we present results for standing 
waves, we will show in § [TO] that the local nonlinear in- 
teraction rates are so fast in both solar-type binaries and 
hot Jupiter systems that the global standing wave as- 
sumption may not be valid; one may then have to work 
instead in the traveling wave limit. 

7.1. Dynamical tide instability for TV < 2 daughters 
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Fig. 6. — Linear energy _Eiin in the modes most resonant with the 
l = k = m = 2 harmonic of the dynamical tide and those modes' 
minimum nonhnear energy threshold _Eth for three-mode coupling 
{N < 2). The results are shown as a function of orbital period 
Paih ior M' = M = Mq. The E^^ points are spaced uniformly in 
Parh'i the density of points indicates the likelihood distribution of 
Eliu near a given Poih- The dotted lines indicate how the various 
quantities scale with Porb- 



For a three-mode system, the dynamical tide is unsta- 
ble if the linear energy of the resonant modes En^ (eq. 
^\ for modes with small Aq) exceeds the parametric 
threshold Eth (eq. [50]) for three- mode coupling to a 
daughter pair. In Figure [6] we show Enn and the mini- 
mum Eth as a function of Porb for M' — M — Mq. To 
find the minimum Eth, we use the analytical expressions 
for K, 7, etc. derived in Appendix |^ and at each Porb 
search the (£, n) parameter space of daughter pairs that 
satisfy momentum conservation and the angular selection 
rules described in Appendix |B] We find that the dynam- 
ical tide in solar binaries is unstable for Porb ^ 40 days 
and that En^ 3> Pth over much of that range. 

The linear energy scales with the mass of the compan- 
ion M' as Ehn oc C/2 oc [M'/{M + M')]^ while Eth is 
independent of M'. The ratio E^n/Eth calculated for so- 
lar binaries in Figure |6] can therefore be used to solve for 
the minimum M' for which the dynamical tide is unstable 
in solar-type stars. We show this minimum companion 
mass M' as a function of orbital period in Figure [71 We 
find that for N < 2 the dynamical tide is unstable for 
Jupiter mass planets M' = Mj out to Porb — 10 days 
and for ^ 10 Earth mass planets out to Porb — few days. 

The analytic scalings shown as dotted lines in Figure 
[6| can be derived using the analytic approximations to 
7, A, t/a and k given in Appendix \X\ The rapid variation 
in Eth and the minimum M' for Porb ^ 10 days are due 
to particular modes coming in and out of resonance. This 
is especially pronounced at short orbital periods since the 
frequency spacing of resonant modes is larger at smaller 

Porb. 

We note that to calculate the minimum M' we used 
the value of Enn given by the lower envelope in Figure 
[6| (eq. [30]). Since this ignores the coincidental possi- 
bility of strong resonances, the minimum M' shown is 
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Period (days) 

Fig. 7. — Minimum mass of the companion M' for which the 
linear dynamical tide is unstable in solar-type stars M = Mq as a 
function of orbital period. We show results for the £ = k = m = 
2 harmonic assuming circular orbits and a non-rotating primary 
(X^™ = 1). The solid line is the minimum companion mass for 
N < 2 while the dashed-dotted line is the min imu m companion 
mass for the collective instability N ^ 2 (eq. I54I V Above the 
dashed line, the local nonlinear interactions are so strong in the 
stellar core that the daughter standing wave approximation is no 
longer valid (eq. 1651 for £i, = 2). The collective instability (A^ 2> 2) 
may invalidate the s tan ding wave approximation for even lower 
companion masses ('i) llO|l . 



somewhat conservative. Furthermore, by approximating 
the dynamical tide with just the most resonant parent, 
we have ignored the driving of daughters by the slightly 
less resonant neighboring parents; we find that including 
these parents decreases the minimum M' by a factor of 
order unity. 

7.2. Dynamical tide instability for N ^ 2 daughters 

In the previous section we assumed that each daugh- 
ter couples to only one other daughter, in addition to 
the parent (either A^ = 1 daughters for a self-coupled 
daughter pair or A^ = 2). However, since \Kabc\ is nearly 
constant for all daughter pairs with \ni, — nd '^ ria (see 
Appendix IB|1. each daughter in fact couples to N ^ ria — 
1000(^Q/2)Pa,io daughters. All the modes are unstable 
if each daughter pair (6, c) £ N satisfies (eq. [35] ) 



27VK„.,Py2 > 



/ 7b 7c - 



A? 

be 



UJhUJc 



(52) 



In order to derive the threshold energy Eth of the system, 
let the radial order n^ characterize the most nonlinearly 
resonant daughter at a given I and consider the N fn ria 
daughters on either side of uq that couple equally well 
to the parent. Of these N modes, the daughter pair 
that maximizes the right hand side of equation (1521) , and 
thereby determines the stability of the system, will be 
the pair that has the largest detuning and damping rate. 
This is the self-coupled pair b = c with nt — uq + na/2] 
it has a detuning |Af,h/a;f,| ~ ria/riQ and a damping rate 
|7h/wb| ~B{no + Ua)^ ^ Bnl, where P ~ 4 x 10-"P6,io 



NONLINEAR TIDES IN CLOSE BINARY SYSTEMS 



15 




rowth 



(N=2) 



M'/M = 10" 



Period (days) 

Fig. 8. — Growth time tgrowthi in units of orbital period, of 
daughters coupled to the linearly resonant parent. The solid 
lines show tgrowth for collectively unstable systems Af 2> 2 with 
M'/M = 2 X lO—* and 10"^. The dashed line shows the max- 
imum tgrowth (multiplied by 0.01) for an Af < 2 system with 
M'/M = 10"=*. 



(see eq. |A14p and the second equality is appropriate if 
^0 ^ -^o = 2. Solving for the no that minimizes the right 
hand side of equation (15^ , we find 



no,: 






1.6 X 10'' 



1/3 



(53) 



with £ = €n * ci (no,*/500)P-\ ~ 32^-^1^. Since 4.* > 



10, f)- 

0,* ± 2, no,* ± Ua) have 



independent of period. This corresponds to a daughter 

^0,* c^ (no.*/500)P-i, 
ia = 2, daughters with (^,n) = 

nearly the same detuning and Kabc as those with {I, n) 
(^o,*,'T-o,* ± no)- Thus a more accurate estimate is iV w 
Zua rather than N Ki Ua- The minimum threshold energy 
is then Etii = (Kafccno,*)~^/24. For a resonant parent at 
the linear energy -Ea.iin of the i — k — m — 2 harmonic 
of the tide, this corresponds to a threshold mass ratio 






> 9 X 10" 



Po, 



1 day 



5/6 



(54) 



where we used equation (|A67p and assumed the primary 
spins slowly relative to the orbit. This minimum mass ra- 
tio for collective parametric driving of daughters is shown 
in Figure [T] Note that this threshold is somewhat con- 
servative in that we again used equation pop for the lin- 
ear energy and thus ignore the coincidental possibility of 
strong resonances. 

At mass ratios just slightly above threshold (by a factor 
of two, say), collectively unstable systems have growth 

1/2 

rates sa 2NKabcEa fl. Their growth rates are therefore 
A^ PS lO'^Pio times larger than N < 2 systems. Thus, 
even though collective systems have a somewhat higher 
threshold (see Fig. [7]), they probably dominate the mode 
dynamics whenever they are unstable. 

In Figure [8] we show the growth times igrowth of the 
fastest growing collectively unstable systems for two com- 



panion masses; for systems just slightly above threshold, 
the daughters all grow on timescales of < 10 orbits. We 
determine tgrowth by numerically solving for the eigenval- 
ues of large networks of daughters (the eigenvalues of ma- 
trix H in § [5]) and using our analytic expressions for ja i 
laim , and Kabc givcu in Appendix |Al To find the daugh- 
ter set with the smallest tgrowth we search over daughter 
i and vary N. For a given i, we couple the N consecu- 
tive daughters on either side of the most resonant daugh- 
ter to each other and the most resonant parent. We in- 
clude daughters with £ ± 2 in the set since their coupling 
also satisfies momentum conservation for ia = ^tido — 2. 
We find that tgrowth decreases with increasing N until 
N w 2na w 2 X 103(4/2)Pa,io. For N > 2na a large 
fraction of the least resonant daughters do not couple to 
each other since they have |An| > n^ and thus Habc — 0. 
We find that these additional daughters are st able. 
Estimating tgrowth using equations (PU]) and (jA67[) . we 

find tgrowth/Porb « 3 X IQ-^ {M /M')P^q^^^ . This is in 
good agreement with the numerical results for mass ra- 
tios sufficiently above threshold. 

8. PARAMETRIC INSTABILITY OF THE 
EQUILIBRIUM TIDE 

In this section we consider the stability of the equi- 
librium tide to parametric driving. Like the dynamical 
tide, the equilibrium tide oscillates at the driving fre- 
quency uj and is thus nonlinearly resonant with short 
wavelength daughters with natural frequencies near lu/2. 
Unlike the dynamical tide, however, the amplitude of 
the equilibrium tide is smallest in the core; at small radii 
Cr.dyn oc r~^ whcrcas ^r,eq oc r^^^ . Since the daughters' 
eigenfunctions peak in the core {^r oc r~^), they couple 
much more weakly to the equilibrium tide than the dy- 
namical tide. We will show that the coupling is in fact so 
weak that the equilibrium tide is only unstable to para- 
metric resonance in solar binaries with Porb ^2 — 5 days. 

A resonant daughter pair (b, c) is unstable to driving 
by the fc-th harmonic of the equilibrium tide if (eq. [35]) 



NT 



where (§ EJ 



n(fc,oq) 

be 



- y^UJb(jJc 



(fc,oq) 
be 
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(55) 
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l+l 



Jcq) 



(56) 



In Appendix [B] we show that the equilibrium tide cou- 
pling coefficient nj^^ ■ (1) is a weak function of daughter 
i and period, (2) has a magnitude of « 0.01 — 0.05 for 
self-coupled modes b — c, and (3) decreases rapidly with 
increasing |An| = jn;, — nd. 

We now estimate the period out to which the equi- 
librium tide is unstable to parametric resonance. Since 
the coupling coefficient peaks strongly for self-coupled 
modes b ~ c, the daughters are not collectively driven 
and A^ = 1. For self-coupled pairs the detuning is 
lAhb/wbl « An^ and the damping rate is "fb/^b — Bn^, 
where numerically we find an average detuning A ~ 0.5 
and P ~ 4 X 10^^^Pio,b (eq. |A14j ). Since the coupling 
coefficient is a weak function of £b and Pb, we can solve 
for the stability threshold by minimizing the right hand 
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side of equation (|55p with respect to ni,. For the repre- 
sentative case of the £ = m = fc — 1 = 2 harmonic in a 
synchronized solar binary (51 ~ 2uji,), we find that the 
equilibrium tide is unstable for orbital periods 



P<1.2(^ days. 



(57) 



There can be brief epochs during which the detuning 
happens to be much smaller than average, i.e. A ^ 
0.5. In the zero detuning limit, the stability criteria is 
determined by the damping rate of the most resonant 
£b = 1 daughter and we find 



P(A6, = 0)<5.l(^^j days. 



(58) 



We have assumed in these estimates that the condition 
|An| = implies self-coupling and thus £i, = £r. In fact, 



pairs with 



f c ± 2 and lAnl — also have k 



(eq) 

be 



0.01 — 0.05 for the £ = 2 harmonic of the tide. In a 
numerical search over potential daughter pairs, we find 
that relaxing this assumption increases the period out to 
which the equilibrium tide is unstable by < 30%. 

9. NONLINEAR INHOMOGENEOUS DRIVING 

There are two mechanisms by which the linear tide can 
drive other modes (see § E]) : parametric driving, which 
we considered in the previous two sections, and non- 
linear inhomogeneous driving, which we consider here. 
Both forms of driving drain energy from the linear tide 
and can thus act as sources of enhanced dissipation rel- 
ative to linear theory. In order to determine the extent 
to which they influence the orbital and rotational evo- 
lution, we need to solve for the nonlinear equilibrium 
(i.e., saturation). In general, the nonlinear equilibrium 
depends on a balance between all three forms of three- 
wave interactions (LLC, LNC, and NNC in §[4l). While 
we defer such a calculation to a future paper, in this 
section we solve for the nonlinear equilibrium in the ab- 
sence of parametric driving and NNC. Such a simplifica- 
tion is justified only if the linearly excited modes remain 
near their linear energies despite their nonlinear coupling 
to other modes. Although we show that this is not al- 
ways the case, this calculation nevertheless demonstrates 
that inhomogeneous driving may be a significant source 
of nonlinear dissipation. 

The steady-state solution to the nonlinear amplitude 
equation in the absence of parametric driving and NNC 
(eq. [20] but ignoring the last two terms on the right 
hand side) is that of a driven oscillator: 



ra(i)=E 



^a [Va + Ka] 



(fe+fc' 



-i{k+k')nt 



■; UJa- {k + k')D, - ija 



, (59) 



where 

[Va + Ka] 

R 



kk 
(k+k') _ f M' 

M 



J2 E Wer,Wi,m'Xt^X'^r' 



tva f m' 



e+i'+2 



^(oq-cq) j^ ^(dyn-dyn) _|_ ^(oq-dyn) 



.(60) 



defined in Appendix IA.6I and their properties described 
in Appendix [Bl This solution (eq. [59]), which is analo- 
gous to that of the linear amplitude equation (eq. [24] ) , 
has contributions from the three distinct forms of driv- 
ing: equilibrium tide-equilibrium tide coupling Ka 



dynamical tide-dynamical tide coupling Ka 



(dyn-dyn) 



and 



equilibrium tide-dynamical tide coupling Ka ■*" . We 
first consider the properties of the solution for modes 
with frequency uJa ^ '^o (high-order g-modes to low-order 
p-modes) and then for modes with frequency Wq 3> wq 
(high-order p-modes). 

9.1. UJa < ^0 

Figure [9] shows the orbit-averaged nonlinear energjij 
Ea^ni = kaP as a function of mode period for each form 
of inhomogeneous driving, as compared to the linear so- 
lution EaMn (eq. [in] and Fig. [T]). We see that for 



Joq-oq) 



Joq-dyn) 



The coefficients Ka , Ka '*'" , and Ka are 



K'a ' " coupling and K^ ^ '^ ' coupling, Ea^nl < EaSm 

over nearly the entire range of modes shown (low-order 
p-modes to n « 10'^ g-modes). This suggests that for 
these modes, the energy dissipation from Ka~ cou- 
pling and Ka ''" coupling is insignificant compared to 
that due to the linear tide. The only exception is the non- 

hnear resonance at Pa ~ Porb/ik + k') due to Ka~ 
coupling. However, since the energy in this resonance is 
below that of the dynamical tide (the linear resonance at 
Pa — Porh/k), it cannot be a significant additional source 
of dissipation. 

For Ka '*'"" coupling, by contrast. Figure [3] shows 
that in the case of solar binaries there are many modes 
for which Ea.ni 3> -Bajin- Furthermore, the energy in the 
nonlinear resonance is '^ 100 times larger than the energy 

in the linear resonancel^ This suggests that Ka ^" ■*" 
coupling may be an important source of nonlinear dissi- 
pation for the dynamical tide, possibly even more impor- 
tant than parametric coupling to daughter pairs. 

It is important to note that Ka '^^ ^" coupling cannot 
lead to enhanced dissipation relative to the mechanical 
power carried by the dynar nical tide (i.e., relative to th e 
upper bound calculated bv lGoodman &: Dicksonl[l998f l. 

This is because Ka '*'"" coupling always occurs below 
the radiative-convective interface (Appendix |B]) and thus 
below the region where the dynamical tide is launched. 
This also implies that if the dynamical tide breaks in the 
core, then the nonlincarly resonant modes cannot actu- 
ally reach the large energies shown in Figure |9l the other 
forms of three-wave interactions (LNC and NNC) must 
become important well before the modes reach these en- 
ergies. 

^ The total energy Ea = |gaP = l^a.lin -I- '"aP contains a cross 
term involving the product of ga.lin ^nd Ta- However, since we as- 
sume in Figure |9] that r^ and Qa.Un oscillate at different harmonics 
of the orbit, the orbit-averaged cross term is zero. In cases where 
Qa lin and ra osciUate at the same frequency (i.e., for particular val- 
ues of k, m etc.), the cross term is constant and the time-averaged 
total energy is not simply |5a linP "I" kaP but can be larger or 
smaller depending on the magnitude and phase of q^ Un and ra ■ 

* Resonances that are coincidentally better than average will 
shift the relative magnitude of the -Ea.nl and _B„ ijn peaks. In Figure 
[9]the linearly and nonlinearly resonant modes both have resonances 
that are typical. 
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Fig. 9. — Orbit-averaged nonlinear energy E^^nl as a function 
of mode period Pa accounting only for the i nhom ogeneous driving 
terms in the nonlinear amplitude equation II20I I. The parameters 
a.re i = I' = m = m' = k - 1 = k' - 1 = 2, M' = M = Mq, 
e = 0.1, and Porb = 10 days. The three contributions to the 



driving — k, 



(eq-eq) (dyn-dyn) 



and 



(eq-dyn) 



— are each shown as 



separate curves. For comparison, the linear energy E^ un is also 
shown (eq. [29]; cf. Fig. [Hi. 

Since i?a, ni ot e^ while -Ea.iin oc e^, we find that in 
the case of hot Jupiter systems, i?a, ni ^ -Ea,iin for all 
modes. We therefore do not expect inhomogeneous driv- 
ing to be an important source of dissipation in planetary 
systems. This difference between planetary and stellar 
companions highlights why it is probably very mislead- 
ing to extrapolate tidal dissipation parameters inferred 
from observations of one set of systems to a very different 
set of systems. 



9.2. 



iOa > Wo 
(cq-dyn) 



1 (dyn— dyn) 

and Ha can 



In § IB. 31 wc found that k 
be significant for high-order p-modes (see the right panel 
of Fig. rTSj) . Since the frequency of these modes is above 
the acoustic cutoff of the solar atmosphere (~ OOwq), 
they do not refiect at the solar surface and cannot there- 
fore form standing waves. Equation (j59p . which assumes 
the driven modes are standing waves, does not therefore 
describe the driving of these high-order p-modes. In- 
stead, one must treat the driving in the local limit; pre- 
sumably the p-modes cannot reach amplitudes anywhere 
near those of equation (|5^ since their group travel time 
across the driving region is much smaller than the driv- 
ing period. We defer the calculation of p-modc driving 
to future work. 

10. LOCAL NONLINEAR INTERACTIONS 

To this point we have assumed that the perturbations 
are all global standing waves. However, a perturbation 
can only be treated as a global standing wave if its growth 
and damping times are longer than its group travel time 
across the star. In this section we show that in nonlin- 
ear theory this is not necessarily the case. The nonlinear 
growth rate of daughters in the core of the star due to 



driving by the dynamical tide can be faster than the rate 
at which the daughters cross the interaction region. We 
specialize to the case of dynamical tide driving since equi- 
librium tide driving, when present, is in the the standing 
wave limit (at least for solar type stars). 

The stability calculations, which are the focus of this 
paper, are independent of the uncertainty of traveling 
versus standing waves that is present in full nonlinear 
theory. This is because the daughter growth rates are 
infinitesimal at threshold, and thus the growth times 
are always longer than the group travel time very near 
threshold. However, the breakdown of the standing wave 
approximation does have implications for the saturation 
of the nonlinear perturbations since it implies that the 
excited waves can grow and damp locally. This, in turn, 
can affect the rate at which energy and angular momen- 
tum are taken out of the orbit. 

Consider an inward propagating traveling wave a 
launched at the radiative-convective boundary re by the 
linear tide. As it approaches the core, this parent wave 
excites daughter waves through three-wave interactions. 
Due to the sharp increase in dni/dlnr in the core (see 
Figs. [T^ and ITU)) . the local nonlinear growth rate Ti,b(r) 
of a self-coupled daughter wave b increases rapidly as the 
parent approaches r ~ . If rbf,(f)ir.f)(^) ^ 1, where 



tr,b[r) 



^bN{r) 



Vr.b{r) 



(61) 



is the time it takes a daughter with group velocity Vr,b ~ 
^b/kr,b to traverse a region of size r and A^ = £b{ib + 1), 
then the standing wave approximation breaks down and 
the daughter must be treated as a traveling wave. In Ap- 
pendix |E] we show that the local nonlinear growth rate of 
a self-coupled daughter b due to three-wave interactions 
with a parent a is 



Tbb(r)c:^2uJbT\qkri 

T 
02 



: 16a;f, 



Eo 



1/2 



A„R 



a-f^^a, 10 



(62) 



where qa is the parent amplitude, T is the angular inte- 
gral given by equation (JA20I) . and the expression in the 
second line is appropriate for r < 0.05R (see eq. |A3| ) . 
The nonlinear interactions are strongest near the par- 
ent's inner turning point r„^in where N(r) — uja — 2a;b; 
the radial group travel time of the daughter across this 
region is tr^bifmhi) ~ '^■^b/^b and therefore at rmin 



Fbbir.bl 



4AbT\qkr^a. 



(63) 



98a;o at r < 0.05R and thus 



In a solar model dN/dr 

rnun/R - 1.3 X lO-^P-,0 

resonant with the linear tide (qa = qa,dyn', eq, 
find that for a synchronous solar binary 



3 X 10 "^P^ In- For the £a = 2 parent most 



Fbbi 



r.6|rn 



2 X lO^AfcePiV^ 



]), we 



(64) 



and the standing wave approximation is invalid unless 
^fcS ^ 10~^. For a slowly rotating star orbited by a 
planet we find 

^bbtr,b\r„i„ ~ 0.3A(, 



M'/M 
10-3 



Pi 



1/6 
10 ■ 



(65) 
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and the standing wave approximation is valid for 
M' /M < 10"^/Af,. The dashed hne in Figure [7] shows 
t his condition for the case h = 2. 

iGoodman fc DicksonI ()1998D computed \qkr^r\a as a 
local measure of the non linearity of waves e xcited by 
the linear tide (see also lOgilvie fc LinI |2007| ) . When 
\qkr^r\a ^ 1, the parent wave overturns the stratifica - 
tion dur ing part of i ts cycle. [Barker fc Qgilvid ()201Clf ) 
(see also lBarkeill2011[ ) carried out numerical simulations 
of gravity waves approaching the core of a solar type 
star. They found th at when the overturn criteria of 
IGoodman fc DicksonI (J1998D is satisfied, the incoming 
parent wave is fully absorbed in the core and therefore 
deposits all of its energy there; if it is not satisfied, the 
parent reflects at its inner turning point and travels back 
out to the stellar surface virtually undamped. Good- 
man fc Dickson's overturn criteria is therefore equivalent 
to a traveling wave condition for the parent. Here we 
have shown that when the parent's amplitude is a fac- 
tor of 4A{,T smaller than this overturning amplitude, the 
daughters it excites are in the traveling wave limit. 

For ^(, < 2 our daughter traveling wave condition is 
very similar to Goodman fc Dickson's parent traveling 
wave condition. For £fc :^ 1, equation (|55|) suggests 
that the daughter standing wave approximation may fail 
even when \qkr^r\a ^ 1, i-e, even when the parent is 
far from overturning the stratification. We also showed 
in § 17.2! that there are collectively unstable daughters 
at mass ratios well below the overturn criteria (Figure 
[7|) with growth rates that are « N/2 ^ 1 times larger 
than that of three wave systems (that of eq. [51]). The 
cases Af, ^ 1 and N ^ 2 both therefore suggest that the 
daughter traveling wave limit may extend to significantly 
lower planetary masses than indicated by the dashed line 
in Figure [7] However, it is not clear whether this means 
that the parent is also in the traveling wave limit, and 
thus fully absorbed, for these lower masses. Assessing 
whether this is the case will require a better understand- 
ing of the daughter driving and saturation in the core of 
solar- type stars. 

We note that both cases (A^ ^ 1 and N ^ 2) involve 
the excitation of a potentially large number of high £ res- 
onant daughters (and since w oc (/n, very high n daugh- 
ters). Capturing these interactions numerically therefore 
requ ires very high spa t ial re solution. This may explain 
why I Barker fc Qgilvid (|2010D do not seem to olaserve ei- 
ther of these effects in their numerical simulations. 

11. SUMMARY AND CONCLUSIONS 

This paper is an initial investigation into the impor- 
tance of nonlinear fluid dynamics for the tidal evolu- 
tion of close binary systems. We derive a formalism 
for including nonlinear interactions in tidal theory and 
describe the physical effects that arise from these non- 
linearities. While the ultimate goal is to understand 
how nonlinear effects may alter the rate of circularization 
and synchronization in binaries, such calculations arc de- 
ferred to a future paper. In this paper, we instead focus 
on determining the conditions under which the standard 
linear theory approximation is (in)valid. We present de- 
tailed results for the tidal forcing of a sun-like star by 
a stellar or planetary companion, but the formalism we 
derive is more general and is applicable to nonlinear tides 
in stars, planets, or compact objects. 



Most previous studies of tides have made the linear ap- 
proximation, accounting for only linear order perturba- 
tions to the background state. Our formalism includes 
the leading order nonlinear corrections. These correc- 
tions have two seemingly different physical origins: (1) 
Internal nonlinear interactions which couple three waves 
to each other. These interactions enable waves that are 
linearly excited by the tide to transfer energy to waves 
that are not linearly excited. This nonlinear coupling 
includes both resonant and non-resonant nonlinear in- 
teractions (§§ HI El fc El). If the linearly excited waves 
systematically lose energy to other waves via nonlinear 
coupling, then the nonlinear interactions can substan- 
tially modify the orbital/rotational evolution relative to 
that predicted by linear theory. (2) External nonlinear 
interactions which resonantly couple two internal waves 
directly to the tide through their gravitational multipole 
moments fii l6.2p . These interactions allow waves that are 
not linearly resonant with the tide to be directly driven 
by the tidal potential. The damping of these nonlinearly 
driven waves can also act as a source of enhanced dissi- 
pation relative to linear theory. 

We have presented two related approaches to solving 
for the nonlinear tidal response of a star. In the first, we 
take the background state of the star to be spherical and 
unperturbed by the companion and we expand the La- 
grangian displacement and velocity associated with the 
tide (both linear and nonlinear) as a sum over eigen- 
modes, with each mode weighed by its amplitudes Qa 
(S I2.1.ip . This extends one of the standard methods for 
studying linear tides to include nonlinear interactions. 
In the second method, we instead take the background 
state of the star to include the linear tidal solution and 
we expand the nonlinear correction to the tidal solution 
as a sum over eigenmodes, with each mode weighed by a 
different mode amplitude r^ fi j2.1.2p . These two meth- 
ods are formally equivalent but we find that each is useful 
for understanding different aspects of tides in stars. The 
second method of solving explicitly for the nonlinear cor- 
rection to the tidal response of a star is particularly use- 
ful for understanding the nonlinear stability of standard 
linear tidal theory (§[5l). 

Since we focus on slowly rotating stars, the stellar 
modes included in our treatment are p-modes and g- 
modes; the latter are particularly important because they 
can have periods comparable to that of the binary sys- 
tem. The orthogonality and completeness of the stellar 
eigenmodes enables us to convert the nonlinear partial 
differential fluid equations into a coupled network of or- 
dinary differential equations describing the evolution of 
each mode amplitude. We now summarize the physics 
contained in these equations, focusing for concreteness 
on the method of expanding the full tidal solution (lin- 
ear and nonlinear) as a sum over eigenmodes. In this 
case, the resulting amplitude equations have linear terms 
and nonlinear terms. Th e linear terms are standard (e.g., 
iPress fc Teukolskvlll977f ) and include driving by the tidal 
force and damping by radiative diffusion^ This is the 



^ While dissipation due to eddy viscosity in convection zones is 
important for the long wavelength equilibrium tide, we ignore it in 
this paper, focusing instead on thermal diffusion damping in the 
radiative zone, which is more important for short wavelength, low 
frequency g-modes in the radiative core. 
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basic physics included in most previous studies of tides, 
with one addition: by solving the mode amplitude equa- 
tions, our formalism allows the modes to be dynamic in 
their interaction with the orbit, rather than assuming 
harmonic response at the forcing frequency. 

The nonlinear terms in the equations of motion are 
parameterized by the coefhcients Kate and Uab (see eqs. 
[To] and [9] in § [2]); these describe the internal (three 
wave) and external (nonlinear tidal driving) nonlinear 
interactions, respectively. Accurately calculating these 
coupling coefficients is nontrivial - these technical details 
are described in Appendix [Aj where we pay considerable 
effort to ensuring that the coupling coefficients can be 
accurately calculated numerically. 

Three wave coupling has been considered exten- 
sively in the l iterat u re, both in stella r seism ology (e.g., 
iDziembowsO Il982t iWu fc GoldreichI |2001[ ) and more 
specifically in the cont ext of tidally excited osci llations 
in binarv svstems ('e.g.. lKumar fc Goodman|[l996D . Non- 
linear tidal driving has n ot, to our knowledge, been s tud- 
ied in any detail before (jPapaloizou fc Pringlelfl98ll con- 
sidered nonlinear tidal driving but neglected three-wave 
coupling which we found cancels strongly with nonlin- 
ear tidal driving). Both types of interactions can lead to 
a variety of physical effects including the non-resonant 
excitation of higher frequency modes and the resonant 
excitation of lower frequency modes (parametric insta- 
bility). We have focused our analysis on determining (1) 
the conditions under which the linear tidal flow is un- 
stable to the parametric instability (§§ [7] & [8|) and (2) 
the efficiency of nonresonant excitation of g-modes and 
p-modes by nonlinear coupling to the linear tidal flow 
(§0. 

In the course of computing nonlinear coupling coeffi- 
cients and growth rates for solar- type stars, we found 
the a priori surprising result that the external nonlinear 
driving of g-modes (Uab) is almost completely canceled 
by a portion of the three wave coupling to the linear 
tide (Kabc)- The implication is that these two effects are 
intimately related. The nature of this relationship can 
be most easily appreciated if we explicitly solve for the 
nonlinear correction to the linear tidal solution (method 
2 in § 12.11) and consider the linear flow as a superpo- 
sition of an "equilibrium tide" and a "dynamical tide" ; 
the former is the nearly hydrostatic part of the linear 
tidal response and the latter is the wave-like, resonant 
part of the response. The three wave coupling coeffi- 
cient that describes the nonlinear driving of a pair of 
'daughter' waves then has a contribution from the equi- 



(dyn) 



librium tide k^^ (eq. [40] ) and the dynamical tide kJ^^ 
(eq. |41|). Since the equilibrium tide contains the vast 
majority of the tidal energy, one might expect it to be 
more pron e to nonlinear ins t ability than the dynamica l 
tide (e.g., iPress et all 119751: iKumar fc Ouataertl [19981 ). 
In fact, we find the opposite: the daughters internal 
driving via three wave coupling to the equilibrium tide 
nearly cancels with their external nonlinear driving by 
the tidal potential (c>c Uab)- As a result, the effective 
coupling of daughters to the equilibrium tide is much 
weaker than their coupling to the dynamical tide, i.e.. 



damentally part of the same process; together they de- 
scribe the nonlinear excitation of daughter modes by the 
nearly hydrostatic response of the star to its companion. 
Another reason for the weak nonlinear driving by the 
equilibrium tide in solar-type stars is that very little of 
the energy of the equilibrium tide resides in the core of 
the star where the low frequency g-modes that the equi- 
librium tide attempts to drive have most of their energy. 
At a technical level, the subtlety of calculating the equi- 
librium tide driving correctly highlights the importance 
of the detailed calculation presented in this paper. For 
example, an order of magnitude estimate of the nonlinear 
driving by the shear of the equilibrium tide {k 



(eq) 
ab 



^ab 



< K 



,(dyn)i 



ab 



(see § 15.21 and Appendices IA.6I and [B ]) . 
Physically, this is because internal nonlinear driving by 
the equilibrium tide and nonlinear tidal driving are fun- 



1) 
overestimates the true driving by al most two orders of 

magnitude (^ 15.31 and Appendix IB. 2 1 Figs. [T3l and fT4]) . 

The equilibrium tide driving of daughter modes, when 
present, is relatively evenly spread throughout the radia- 
tive zone (left panel of Fig. [T3]); by contrast, the driv- 
ing by the dynamical tide is highly concentrated in the 
core of the star (Fig. [T2]) . Hence the two mechanisms 
for seeding fluid instability are physically quite different. 
Quantitatively, we flnd that the equilibrium tide is only 
unstable to nonlinear driving of resonant g-modes in so- 
lar binaries if the orbital period is < 2 — 5 days (§ [5]). 
This implies that instability of the equilibrium tide in 
slowly rotating solar-type stars cannot help explain the 
observed circularization of solar binaries out to orbital 
periods of ~ 10 — 15 days. An interesting remaining pos- 
sibility is that the equilibrium tide in rotating solar-type 
stars may efficiently couple nonlinearly to inertial waves. 
In particular, inertial waves in solar type stars have most 
of their energy in the outer convection zone where much 
of the energy in the equilibrium tide also resides. 

Recently, the elliptical instability has been invoked 
as a source of nonlinear driving in close bin ary sys- 
tems (iLe Bars et all 120101: iCebron et al.|[201lt see also 
IPress et al.lll975l:lSeguinlll976D . The elliptical instability 
is clearly related to the equilibrium tide instability con- 
sidered here; both instabilities involve the destabilization 
of internal modes of oscillation by the long wavelength, 
hydrostatic, tidal perturbation. They are not equivalent, 
however, as the elliptical instability drives inertial modes 
whereas we focus on g-modes. 

For solar type stars, we find that the dynamical tide 

driving rates of short wavelength g-modes (ex Hab^'^ ) ^^^ 
orders of magnitude larger than the equilibrium tide driv- 
ing rates. In order to determine when linear theory is 
valid, we computed the parametric instability growth 
rates over a range of orbital periods and companion 
masses (§[7l). We find that for orbital periods from days 
to weeks, the dynamical tide part of the linear tidal flow 
is unstable, even for companion masses much smaller 
than a Jupiter mass (see Fig. [7]). The linear tidal so- 
lution in the star is invalid even for a 10 Earth mass 
planet if the orbital period is < 5 days! The degree to 
which the true tidal flow differs from the linear one de- 
pends on the saturation of these nonlinear instabilities, 
which will be investigated in a future paper. This ques- 
tion is of particular current interest given the large num- 
ber of low mass planets orbiting solar-type stars being 
discove red by transit surve ys such as the Kepler tele- 
scope (jBorucki et al.l [20111 ). The orbital properties of 
these systems are likely shaped in part by the nonlinear 
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instabilities discussed in this paper. 

In our study of the parametric instabihty of the dy- 
namical tide we found that it is subject to a "collective" 
version of the parametric instability. In the literature, 
parametric growth rates have traditionally been derived 
by considering how a single pair (N < 2) of coupled short 
wavelength daughters interact with the time-dependen t 
background of a parent mode (e.g., iDziembowskil 1 1 982l ) . 
We have extended this analysis to allow for A^ ^ 2 
daughters, each of which is coupled to the parent and the 
other daughters. We find that for solar type stars, large 
sets of daughters are coherently driven by the dynamical 
tide {N w 10^ — 10'^ for periods P ^ 10 days and com- 
panions M'/M > few X 10^^). When collectively un- 
stable, the daughters all grow as a single coherent unit, 
i.e., maintaining phase coherence. This coherence is self- 
consistently generated by the nonlinear interactions even 
if it is not present in the initial amplitudes/phases of the 
modes (Fig. U). In the collective parametric instabil- 
ity, the growth rates of the daughters are significantly 
larger than in standard three wave coupling, by a factor 
of « N/2 (Fig. [5]). As a result, the collective instability 
may have important implications for the nonlinear damp- 
ing of the dynamical tide. More generally, it will also be 
important to revisit previous astrophysical applications 
of parametric instability to understand the consequences 
of the restrictive focus on three-mode coupling, which 
does not necessarily capture the parametric instability 
threshold or fastest growth rate. 

In th eir study of the dynamical t ide in solar-type bi- 
naries, IGoodman fc DicksonI (jl998[) noted that short- 
wavelength gravity waves traveling inwards from the 
radiative-convective boundary (where they are excited) 
reach sufficiently large amplitudes that the waves 
"break " at the ce n ter o f the star. [Barker fc Qgilvii 
(|2010| ) and iBarkeii ()20lH ) confirmed this result numer- 
ically. For hot Jupiter companions with masses < 
2Mjupitor, however, they found that the incoming grav- 
ity waves d o not break in the core of the star (see also 
IQgJlvie fc Lin 2007). 

In an effort to understand what happens in the regime 
where the incom i ng wa ve does not break in the core, 
[Barker fc Qgilvid (J2011I ) performed a stability analysis 
of a non-linear standing internal gravity wave in the cen- 
tral regions of a cylindrical "star" . Although their meth- 
ods differ from ours, they obtain (non-collective) daugh- 
ter driving rates simila r to ours (see § 15. 4p . Like us, 
[Barker fc Qgilvid (|2011[ ) do not attempt to solve for the 
saturation of the unstable daughters and thus do not 
explicitly solve for the tidal Q of the star. However, 
they argue that the driving rate of the fastest growing 
daughter mode places a lower bound of Q > 10^. We 
are not convinced that this result is correct for two rea- 
sons. First, both our study and theirs find that the par- 
ent excites many daughter pairs, many of which have 
(non-collective) growth rates just slightly below that of 
the fastest growing mode. These excited daughters will 
each drain energy from the parent. Since there are so 
many such rapidly growing daughters, it seems possi- 
ble that the parent will lose more energy to all of these 
modes th an it does to the s i ngle f astest growing mode. 
Secondly, [Barker fc Qgilvid ()2011f ) do not seem to ob- 
serve the collective instability, and thus might be un- 



derestimating the rate at which daughters drain energy 
from the parent F^ Qn e possible explanation for why 
[Barker fc Qgilvid ()201lD do not see the collective insta- 
bility is that their diffusivity, which they chose in order 
to obtain numerically converged solutions, is set so high 
that it artificially suppresses collective driving. In §[7] we 
show that the companion (planet) mass above which the 
collective instability sets in lies well above the thresh- 
old for the standard three-mode instability (but still well 
below the wave breaking threshold; see Fig. [7]). Thus, 
if their diffusivity is set too high compared to the true 
value in a solar-type star, they might be in the regime 
where standard driving is unstable but collective driving 
is artificially stabilized. 

In addition to the driving of resonant high i daughter 
waves, there is efficient non-resonant three mode coupling 
(nonlinear inhomogeneous driving) between the dynami- 
cal tide, the equilibrium tide, and a wide range of stellar 
p- modes and g- modes (§[2] and Figs. [n]&[T5]). How this 
coupling saturates is, however, unclear. Even if these 
alternative nonlinear couplings turn out to be more im- 
portant than wave breaking in damping the dynamical 
tide, they are unlikely to modify th e orbital evolution due 
to the dynamical tide calculated bv [Goodman fc Dicksorj 
(|1998f ): this is because the linear energy input rate into 
the dynamical tide is independent of its damping rate in 
the limit of efficient damping. Nonetheless, identifying 
the correct damping mechanism for the dynamical tide in 
solar binaries is important because it determines where 
in the star the tidal energy is dissipated. 

In some cases, we find that the parametric growth 
times due to driving by the dynamical tide can be so 
short that the excited daughter waves do not have time 
to travel through the core of the star (at the group ve- 
locity) before their amplitudes are expected to become 
highly nonlinear (^ llOp . This is particularly true for the 
"collective" parametric instability (Fig. 10). In this case, 
the standing wave approximation employed here may be- 
come inaccurate and a traveling wave point of view may 
be more appropriate. While this does not alter our con- 
clusions about the stability of the linear solution, it will 
likely affect how the nonlinearities saturate and thus the 
rate of tidal dissipation. 

To conclude, we note that the general formalism devel- 
oped in this paper can be extended to the tidal forcing 
of other types of stars. For solar-type stars, the reso- 
nant g-mode amplitudes are large at the center of the 
star, where the equilibrium tide amplitude is small. In 
other stars, resonant wave amplitudes may be large at the 
stellar surface, where the equilibrium tide is also large. 
Equilibrium tide driving may thus play a more prominent 
role in those stars. 



The authors would thank Peter Goldreich, Pawan Ku- 
mar, and Yanqin Wu for useful conversations during the 
development of this work and the referee for detailed and 
helpful comments that improved the manuscript. This 
work was supported by NSF AST-0908873 and NASA 

^^ In the abs ence of collective driving, the driving rate Ffj^ oc 
ui/n (see § 15.41 1. Collective driving increases the driving rate of 
each daughter by a factor of approx imately N Ki n and ther efore 
the collective driving rate F;,;, oc oj. [Barker fc OgilvI3 120111 ) find 
Tiji, oc uj/n even above the collectively instability threshold. 



NONLINEAR TIDES IN CLOSE BINARY SYSTEMS 21 

NNX09AF98G. P. A. is an Alfred P. Sloan Fellow and was supported in part by the Miller Institute for Basic 

received support from the Fund for Excellence in Science Research in Science, University of California Berkeley, 

and Technology from the University of Virginia. E. Q. and the David and Lucile Packard Foundation. 

APPENDIX 

A. COEFFICIENTS IN THE AMPLITUDE EQUATION 

In this Appendix we derive expressions for the coefficients of the two forms of amplitude equation (eqs. [6] and |20|). 
In § lA.ll we describe our calculation of the linear tidal displacement ^ y,„ ne eded to calculate some of the co effici ents. 
Expressions for the linear-order coefficients 7a a nd U g are given in § IA.2I for the nonlinear tide Uab in § IA.31 and 
for the three-mode coupling coefficient Kabc in § IA.4I Analytic approximations to Kgb c are given in § IA.4.11 and the 
correction t o Kghr. when one of th e modes is ^ii„ rather than an eigenmode is given in § IA.5I 

Following 'Schen k et al.l (|2002| ). these coefficients are found by expanding the displacement ^ as a sum of linear, 
adiabatic eigenmodes (eq. [5]) and contracting each mode ^g with the linear and nonlinear internal (pressure, buoyancy, 
perturbed gravity) forces fi and /2 and external (tidal) forces atido (eq. [3]) according to 

{^^,F)^l-Jd^xpC-F, (Al) 

where F is one of the forces f^, /j, or patidc and we choose a normalization in which each mode has the same energy 
Eq = GM^/R. The normalization integral for mode a is then (eq. [7]) 

Eq = 2x (fx( -ujIp 11^ I^ + potential energy | (A2) 

= 2X 2x1^2 1^3^^ |^J2 

c:,2ulAl fdrr'p{r)eg.ir), 



where one factor of two arises since both a mode and its physically indistinct complex conjugate mode are included 
in the energy (see eq. [5] ) and a second factor of two arises from the fact that the potential energy equals the kinetic 
energy. In the last equality we specialized to the case of low frequency g-modes for which the horizontal displacement 
^h is much larger than the vertical displacement £,r- This approach to solving the nonlinear partial differential equation 
(j4t is sometimes referred to as the Galerkin method, a particular implementation of the method of weighted residuals 
in which the weighting functions are chosen to be the eigenfu nctions of the linear system ( s ee, e.g., [Finlavson 1972). 
The expression for the linear force fi[^] is standard (see, e.g.. lLvnden-Bell fc Ostrikeiill967l : iSchenk et al.ii2002i ). The 
linear eigenfunctions ^g and frequenc ies ujg are found by solving t he eigenvalue problem /i[^a] = —pt^a^a using the 
Aarhus adiabatic oscillation package (jChristensen-Dalsgaardl [20081 ) . We use a 5 Gyr old solar model taken from the 
EZ code ([Paxton [2004) . The express ion for the leading order nonlinear interaction forces /2K,^] and the nonlinear 



Schenk et al.l (|2002D and given by the nonlinear terms in their equations (4.8) and (4.9). 



tide is derived in 

In Figure [TOj we show the radial profile of the g-mode eigenfunction that is linearly resonant with the tide for 
Porb = 10 days for a solar model. Also shown is the WKB approximation to this high-order mode (see, e.g.. lUnno et al.l 

Cr{r)=A{r) am (p{r), (A3) 

a(r)-A(r)-^cos^(r), (A4) 

where A = £{£ + 1), kr = AN/cor is the radial wavenumbcr of the mode, (j){r) ~ J^ krdr + tt/A is its phase, 

and AP = 27rV/o''' Ndlnr ~ 1.27/wo- For r < 0.05i?, N/r ~ 98ujo/R and p ~ 150 g cm'^ ~ 26 M/R^, and to a 
good approximation A{r) ~ 3.5 x 10*cm(r/i?)^^ ~ (_R/200)(r/i?)^^. This WKB expression is valid in the propagating 
region (i.e., where uj < N) and, as we show below, is useful for estimating the magnitude of some of the coefficients in 
the amplitude equation. 

A.l. Linear tidal displacement ^j;,j 

In order to calculate some of the coupling coefficients (e.g., n^^ and Hf,/'^ defined in eqs. [40] and [41]) we must 
calculate the equilibrium tide and dynamical tide displacements ^gq and ^^yn- We do so by numerically solving the 



22 



WEINBERG, ARRAS, QUATAERT, & BURKART 







T3 



O 




R 



R 



Fig. 10. — Left panel: Radial eigenfunction ^r {dark gray line) and horizontal eigenfunction ^i^ (light gray line) of the mode {ia^na) = 
(2,486) for a solar model. This mode is linearly resonant with the tide for Porb = 10 days (see Figure [Jl. The solid black lines show the 
WKB approximation to the mode (eqs. [A3 1 IA4| without the sinusoidal factors). Right panel: Radial profile of the three- mode coupling 
coefficient dKahc/dlnr {solid black line; divided by the angular in tegral T) of a parent with {£a,na ) = (2 , 486) and a self-coupled daughter 
b = c with (<?(,, ni,) = (2, 971). The full Kabc expression IIA55IA62I I is dominated by the seco nd te rm liA56l l. plotted here as grey disconnected 
points. The dashed line is the WKB approximation to the amplitude of dK^i,^/dlnr (eq. IA64n . 

linearized momentuni, mass, and energy equations governing the adiabatic response to a perturbing potential U 

piiin = -V6p + gSp-pV{d(P + U) (A6) 



Sp=-V ■ (p|ii„ 
pN^ 



Sp 
op^^ 



9 



ir, 



lin- 



(A7) 
(A8) 



Here ^ij„ = ^on+€dyn is the Lagrangian displacement vector of the linear tide, 5 indicates an Eulerian perturbation, 
and the other symbols have their usual meaning. We can separate the radial, angular, and time dependence of the 
fluid variables in the usual way. Defining the potential ip = 5p/ p + 5<f> + U , the radial momentum, continuity, and 
Poisson equation give 



dr 

'^^r.lin 

dr 



N^ 



{il; - S(f> - U) - {N^ 



A?i! 



^'Z e.,1 



- ^ )^rAu 
Ip — S(j> - 



u 



dr r r^ 



A-kGp 



i;-S<f)-U 



9 



(A9) 
(AlO) 
(All) 



where 



d6(j)/dr. Given i/;, the horizontal displacement is S,h, 



lin = 'pjruP'. At the center we impose the regularity 
conditions Cr.iin = ^Cft.,iin and (5(/) oc r^. At the surface, the condition d54>/dr = —{i+ \)54)/r picks out the solution that 
decreases outward and we require the fluid to be hydrostatic by imposing 6p = pgCr.iin- 
We define the equilibrium tide displacement ^ as the solution to equations (JA6IIA8P in the limit lo -^ Q. We fincF^ 

^r,oq = ~ : (A12j 



Spcq = P9^r,cq, aud 5 Pcq = -^r,cqdp/dr; wc computc 
solving for ^ii„ and l^q we find ^dy„ = ^ii„ - ^^^. 



9 

1 d 



{rX: 



cq; 



(A13) 



rA^ dr 

oq by numerically solving Poisson's equation (JA11[) . Upon 



^^ These expressions are valid as long as A'^^ 7^ 0; they therefore 
apply in radiative regions {N-^ > 0) and convective regions {N'^ < 



0). Given our definition of ^ , we have ^^y^ 7^ in a convective 
region (see lGoodman fc Dicksonlll998l : ITerauem et al.lll998) . 
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Fig. 11. — Left panel: The top plot shows the £ = 2 Hnear overlap integral Ia2m as a function of mode frequency in units of uiq = 
(GM/R^)^''^ for a solar model. The solid line is the numerical result, while the dashed line is the analytic approximation (eq. lAlGI ). The 
bottom plot shows the damping rate due to radiative diffusi on. Li nes representing £ = 2 — 8 modes are shown; they very nearly lie on top 
of each other after the A^ scaling has been removed (sec cq. IA14n . Right panel: Radial profile of lalm for b. solar model. The parameters 
are £ = 2 and {£a, n-a) = (2, 485), which corresponds to a mode period of Pa ~ 5 days. 




A. 2. Linear damping rate ja o-nd linear tide coefficient Ua 

The temperature fluctuations that accompany the density perturbations of the excited g-modes are smoothed out by 
radiative diffusion. This is the dominant linear damping mechanism for high-order g-modes. We calculate its rate jg 
by co mputing the quasi-adiabatic work integral (see, e.g.. lUnno et al.|[l989HTerquem et al]|1998t [Goodman fc Dicksoiil 
Il998f ). The results, shown in the bottom left panel of Figure [TTl are well- fit by the formula 



7„(^,c^)~2x10-"A2 



UJQ (X k 

w, 



2 
r,a' 



(A14) 



Since the radiative diffusion is proportional to the second derivative of the temperature fluctuation, the damping rate 
of short-wavelength perturbations scales as fc^ . 
The overlap integral 



1 



MR<^ 



J d^xpC ■ V {r%„,) ^J^J drr'^Hpg 



2^+1 (50(i?) 
47r GA//i?' 



(A15) 



upon which the linear tide coe fficien t Ua depends, is shown in the top left panel of Figure [TT] for i — 2. The second and 
third expressions in equation (jA15|) can be derived by integrating by parts (see lZahnlll97Cl[ ). The radial profile of I aim 
for a high-order mode is shown in the right panel of Figure 1111 Because the linear overlap laim involves an integral 
over a single high-order mode, it is highly oscillatory from r = out to the radiative-convective interface at r ~ O.Ti?©; 
the oscillations cancel almost perfectly throughout this region. The main contribution to laim comes from a region 
near the radiative-convective interface; since the mode becomes evanescent in this regi on, its wavelength becomes very 
long, allowing it to finally couple well with the large scale tidal potential (|Zahnlll9"75l ). 
For mode periods Pa ^ 1 day. 



^alra 



2.5 X 10-3 



Wo 



11/6 



(A16) 



I Goodman fc DicksonI (|1998[ ) also found an u?-^!^ scaling by matching approximate solutions of the inhomogeneous 
linear fluid equation across the radiative-convective boundary. As a check of the numerical accuracy of our calculation, 
we verified that our overlap integrals satisfy the sum rule 



J2^lla2m^5G dxpx\ (A17) 

a '' 

to high accuracy, where x = r/R (see lReiseneggeiill994l : note that he uses a different eigenfunction normalization). 
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A. 3. Nonlinear tide coefficient Uab 

We now calculate the nonlinear overlap integral JaUm (eq- |13]). upon which the nonlinear tide coefficient Uab 
depends. Using the covariant basis with vectors e^ — hid where hr — l,h0 — r,h^ = rsm6, the components of the 
Lagrangian displacement vector are 

ah dYa ah dYa ' 



where a^ = ar{r) and a^ — ah{r), and Ya = yi^,ma{^ 
then consists of the following terms: 



dj^i d 



(A18) 



r de 'rsin^e dc, 
is the spherical harmonic function. The integrand of Jabim 



Caa{r'Yt).,rr=i{i' 



ly-'OrbrYaYbYt 

00 00 

1 ^ 5n ay* 

sin^ e ' 






Ceb{r'Y,).,oe = r'-Wbh 



OYg OYb 

00 00 



O'Yt 
00^ 



eaCyY,),o^ = r'-'ahbh 
^UtyY,),^^ = r'-^ahbh 



1 OYg OYb 

sm^0 00 Ocj) 

1 OYg OYb 

sin^0 04) 0(f) 



iY 



O^Yt cos 0OYt 



sm( 



T5- +sm0cos^— — 



a^r* 



'sin2 6'r, 



(A19) 



where £ and Yt refer to the harmonic of the tidal potential and the 0r, (pr, and 40 terms not shown are analogous 
to their symmetric count erpart. The angular integrations can be done analytically. Following the notation of Wu & 
Goldreich (2001; see also lKumar fc G oodman , lOgG. ). we define 
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dflYr, 
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*-a ^b ^c 

nia rub rUc 



1 OYb OYc 



T 



(A? 
T 



K-a ^h ^c 



A^-A^l 



dn [5™5''"(v.v,y„)(v™n)(v„y.)] = ^ [a^ - (a^ - a^) 



(A20) 

(A21) 
(A22) 



where <?'■' is the metric tensor on the unit sphere, the matrices on the right hand side of equation (|A20p are Wigner 
3-j symbols, and other terms (e.g., Fb, Fc) are found by permuting indices; note that these expressions apply to any 3 
spherical harmonics denoted by a,b, and c. For the particular case of calculating Jabirm "c" here stands for the tidal 
potential labeled by £ and m. With these definitions we find 



<^ai 



1 



bhn - 



1 



dr pr^ [£{l - l)a^6^r + {£ - l)arbhFa + {£ - l)ahbrFb 
+ahbh {Gt + iFt)] . 



(A23) 

The angular integral T is subject to the selection rules \£b—£c\ < £a < £-b+£c with £a+£b+£c even and ma+mb+nic = 0. 
These rules, which ensure conservation of angular momentum during the nonlinear interactions, allow modes with 
(■a^^b 7^ 2 to couple to the £ = 2 component of the tide. The linear overlap integral lainn by contrast, involves a 
product of two spherical harmonics and thus vanishes unless the angular degree of th e mod e equals that of the tide. 
In the propagating zone of low frequency g-modes, K^h 3> fr and the last term in (|A23[) is dominant. However, as 



bc£m' 



we explain in § 15.21 this term cancels with the inhomogeneous piece of the three- wave coupling coefficient n 

A. 4. Nonlinear coupling coefficient Kabc 

In this appendix we derive the form o f the thre e-wave coupling coefficient. We follow the treatment of Wu & 
Goldreich (2001, hereafter WGOl; see also iWul 1 19981) . who found that significant cancellations between large terms in 
the coupling coefficient can lead to decreased numerical accuracy. They showed that many of the cancellations can 
be removed analytically using integration by parts, resulting in a mathematically equivalent expression that can be 
integrated with much higher numerical accuracy. 



NONLINEAR TIDES IN CLOSE BINARY SYSTEMS 



25 



The final result we derive below differs in four respects from WGOl. First, their equation (A15) should be multiplied 
by —3/Eq to obtain our normalization. Second, their equation (Al) neglects a term involving the gravitational potential 
of the unperturbed background (this term is not the perturbed gravitational potential that one can justifiably ignore 
when the Cowling approximation is valid). We find that this term is large, and leads to some cancellation in their 
equation (A15). Third, we believe there are errors in terms 5, 7 and 9 of their equation (A15), as discussed below. 
These errors are also contained in their equation (A14). Fourth, we include the perturbed gravitational potential, 
both in the original form of the 3-wave Lagrangian and when using the linear equations of motion to simplify the 
expressions. Lastly, for clarity, we write out all permutations of mode indices explicitly. 

O ur fina l expression for Kabc is given by lines (|A55IIA62|) . While quite complicated in general, we show in ?? IA.4.1I 
and lA.4."2l that the expressions simplify greatly when the daughters are short wavelength modes. In particular, we show 
that dKabc/dhir approximately equals the energy density of the daughters dE/dr times the shear of the parent kr^r- 
This is true regardless of whether the parent is a short or long wavelength wave (e.g., dynamical tide or equilibrium 
tide). 

For brevity, in this subsection we use the notation that the displacement vector for a mode a is a, with radial and 
horizontal components Ur and ah, an d similarly for modes b and c. The three-wave coupling coefficient is given by 
equation (4.20) of .Schenk et all (|2002[ ): 



Kabc=7r7T- I d^x p 

ZhQ 



(ri-i)2 + 



d\np- 



V • aV • 6V • c 



-(Fi - 1) (a^6^ V • c + 6^c^ V • a + c>^, V • b) + c^^b{j!^ + al^c'.b'; 
-—— / d^xp [a'VS(t)c-aj + Vc'5(j)a:ij + da^54)b.^j + a*6^c''0;yfe] • 



(A24) 



This expression is symmetric with respect to the three modes. The semicolon denotes a covariant derivative. 

To evaluate this expression, WGOl first performed the angular integra tions. Using WG Ol's notation, these integra- 
tions can be expressed in terms of the angular integrals T, F^, Ga (eqs. |A201IA2TllA22| '). and 



Va = AlAlT-Fa-S. 



KFc) 



(A25) 
(A26) 



The selection rules on the angular integral T (§ lA.Sp enforce angular momentum conservation during the 3-wave 
interactions. 
After performing the angular integrals, WGOl showed that many terms cancel, and several terms could be integrated 
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by parts. Our expression for the remaining terms, using WGOl equation (A5-A8), is 
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The additional terms not included in WGOl are the ''a^b'c''(t)-ijk" terms on line IA30I arising from the background 
gravity, and the "a^VS4>c;i,j terms on line IA3l1 arising from perturbed gravity. 

To this point, the linear equations of motion have not been used. We now follow WGOl and integrate by parts, 
simplifying the resulting expressions using the equations of motion. While WGOl's equation (A9) still holds, we must 
include the perturbed gravity as follows: 



d^r 2 ^ A2 
ar r r 

TipV ■ ^ = pgir - oj^pri 



6. + V • I 



d , „ .^ A2 ^ 

— (TipV • $) = —pgih 
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2 , 25 dg\ . 

r dr I ar 



(A32) 
(A33) 
(A34) 



Note that equations (jA33[) and (IA34P are the homogeneous equations of motion (i.e., they do not include inhomogeneous 
terms involving the tidal potential U) and are therefore only appropriate if all three modes are solutions of the 
homogeneous equations. As we show in § IA.51 when one of the modes is replaced by a solution of the inhomogeneous 
equations of motion (i.e., the linear tidal displacement ^ii„), there are additional terms in the final expression for the 
coupling coefficient. 
The largest terms on line IA29] have the form "(i/(ir(aft,6/j)". Integrating these three terms by parts, plugging in 
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equation (IA34p . and canceling terms against the last term ''—2SpgahbhCh" on line I A28l we find 
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(A35) 



where we ha ve us ed equation (jA25p in the cancellation. 

Next, line IA291 contains terms involving d^h/dr. These terms can be integrated by parts and then simplified using 
equation (IA34[) . The derivatives of the form d^r/dr can then be eliminated from line IA29I using equation (jA32[) . The 
resulting expressions contain derivatives only in the divergences V • ^. For terms containing just one factor of V • ^, 
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we use equation (|A33I) . After significant cancellation of terms, and collection of like terms, we find 
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(A36) 
(A37) 
(A38) 

(A39) 
(A40) 
(A41) 

(A42) 

(A43) 
(A44) 



(A45) 



(A46) 



We find three differences between our expression and WGOI's equation (A15). First we find that their term 5 
should be — 3ripV • ab^CrT, instead of — (SFi + f)pV • abrCrT. Next, their term 7 should be —pgarbrCrT/3 instead 
of 2pr~'^arbrCrT. In their term 9, we find 5pg — 2p/r should be replaced by 4:pg. Lastly, inclusion of the background 
gravity terms found in line I A44l which are not contained in WGOI's equation (Al), cancel line IA38I Hence neglect of 
the background gravity terms is not justified. 

Partial cancellation occurs between t he d4 > terms arising from integration by parts (line IA45P and the terms found 
in our equation (jA24[) contained in line IA461 In particular, the large terms of the form pahbfiS(j)cGc cancel out. 

The coupling coefficient as written on lines IA36IIA461 significantly reduces the cancellation error compared to the 
original form (our eq. |A24| ) . Yet near the center of solar- type stars we find that there is still considerable cancellation 
between the terms of type a^ferV • c, arb^Cr, and ahbrCr- These terms rise the fastest toward the center, and come 
to dominate there. However, after integration they cancel against each other. For high radial order g-modes, these 
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individual terms can be orders of magnitude larger than their sum. We transform these terms as follows: 
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(A47) 



(A48) 



(A49) 



(A50) 



Going from line IA47I to line IA48[ we plugged equatio n (|A34[) into the first three terms and equation (|A32I) into the 
last three terms. This expression was simplified in line IA49I Going from line IA49I to line lASOl we integrated the term 
d/dr{arbrCr) by parts and simplified. The integrand is now not so steeply peaked toward r = 0. Further, the last 
terms of the form ahbrCr are now proportional to w^, not the large frequency g/r. 

The terms involving 6(1) can also be further simplified in order to remove terms that cancel against each other near 
the center. The Sdr terms can be transformed as follows 
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(A51) 

(A52) 

T (A53) 

(A54) 



To get the first equality we used (JA21I) . to get the second we used equation (jA32p . and to get the last we integrated 
by parts. The 5(f>a and S(j)b terms yield analogous expressions. 
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^abc 



dr 



Tr^p ri(ri + i) 



We thus obtain our final result: 
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Since dp/ dr can be written in terms of iV, no numerical derivatives are needed in this expression. 

In the right panel of Figure [TU] (see also Fig. [T^ we plot dKabc/d^^f as a function of radius for a parent a resonant 
with the linear tide and a self-coupled daughter b = c with period Pb ~ 2Pa = 10 days. 

A. 4.1. Analytic estimate of Hahc for three high-order g-modes m solar-type stars 

We now derive an analytic estimate of Kabc for the coupling of three high-order g-modes (see a lso 'DziembowskH 
119821) . Over such a mode's propagation region (i.e., where w < N), the WKB approximation (eqs. |A3[ IA4| ) is very 
accurate (see left panel of Figure FlOl) . Over most of this region uj <^ N/A and thus ^/j ^ ^r- Below the inner turning 
point in the core (w > N), the mode can be approximated with Bessel functions. The gravitational perturbation 6(j) 
due to the mode is negligible throughout the star (i.e., the Cowling approximation is appropriate). An inspection of 
our f inal Kabc expression (lines IA55tlA62]) reveals that as a result of these properties, the second group of terms (line 
IA56[) is the largest of the 8 groups of terms by at least a factor of a few throughout all but the innermost part of the 
star (see grey points in right panel of Figure llOp . 

First consider the region where to <C N/A for all three modes. Using the WKB approximation (eqs. |A3[ IA4) ) we 
find 
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(A63) 



(A64) 



Here H is the pressure scale height and we used equation (jA33p . The brCh and b^Cr terms in equation (jA56[) are 
negligible since they involve products like sin(/ifccos0c — sin(0(, — 4'c)/2 and thus very nearly vanish (the daughters 
must be similar by momentum conservation and thus (pb — 4>c)- The brCr term, by contrast, involves the product 
sin (j!)f, sin 0c — cos((/)b — 4'c)/2. The factor g/TiHN'^ varies from ~ 1 — 4 between the center and convective zone of 
a solar-type star. The expression on the second line gives only the component that varies slowly with background 
quantities and not the high-frequency oscillatory component cos0a sini/);, sin(/)c- As shown by the dashed line in the 
right panel of Figure 1101 equation (JA64I) is a very good approximation to the amplitude of dn/d In r over nearly the 
entire star. Note that since the local energy of a mode within a shell of thickness dr is dE ~ pr^N^^^dr we see by 
(JA63P that insofar as br ~ Cr, 

dKgbc _ djEb/Ep) dar d{Eb/Eo) 

1 :— ~ 1 Kralr, (AOOJ 



dlnr 



d In r dr 
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i.e., dn/dhir is proportional to the parent's local shear dar/dr; we will use this expression to evaluate the local 
three- wave coupling in Appendix lEl 

Now consider the region in the core between the parent's inner turning point r^ ~ uJa{dN/dr)^^ and the daughters' 
at rfi^c — ''a/2. The daughters are again given by the WKB approximation while the parent is evanescent and given 
by Bessel functions a^ ~ Axji^{x) and Uh — {rA'^)~^d{r'^ar)/dr, where A is given by equation (jASP and x — kr^af- 
For X ^ 1 and ta = 2, J2 — x'^ /8. Plugging in values corresponding to the core of a solar model [dN/dr ~ 98ajo/^j 
Ta ~ 1.3 X 10-4p-ijji?, p ~ 150 g cm-3) we find 
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d In r d In r 
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term— 2 
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In the standing wave limit (see ? llOp . the global coupling is dominated by the region ?'q/2 < r < r^ and for resonant 
coupling Pa ~ Porb/2 we find 



in good agreement with the full numerical calculation. 

A. 4. 2. Analytic estimate of Habc for the equilibrium tide coupled to two high-order g-modes in solar-type stars 

We now consider the case where the parent a is the equilibrium tide, which corresponds to the case Kabc = i-^bc 
(see §§ 15.21 and lA.Sp . We then have Or ~ {r/RY'^'^R and h^an « (^ + A)ar at r ^ R a nd on e c an sho w that six of the 
eight terms in Kabc are of similar magnitude (the exceptions are terms 1 and 7 [lines IA55I and IA61| , which are much 
smaller and c an be ignored). For example, since V • a = and V • 6 ~ br/TiH for short wavelength daughter g-modes, 
term 2 (|A56|) is approximately 

Tg pr^N%rCr. ^^ Tg£ djEb/Ep) f Cr\ ^ ,,__. 

-^^r ~ ^,^ ,,,,, :r ^ Aa,. (A68) 
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2TiHN^ Eo 2ATiHN^ dr \br 



The ratio Cr/br implies that kj^^ is small unless bj. ~ c^ (i.e., \nb — nd ^ 1), since otherwise the integrand oscillates 
rapidly about zero on the scale of the long wavelength parent. The other five terms are of similar magnitude (although 
their signs can vary). Since g£/2ATiHN^ w 1 throughout the radiative zone where the daughters propagate, we find 

dj^ ^ M/Eo)^^^ ^ MEo)^ ry-^ (^gg^ 

dmr dr dmr \R/ 

where the second equality applies at r ~ i?. Therefore, as we found in the case of a high-order parent (eq. jA65j ). the 
local coupling of daughters to the equilibrium tide is proportional to the local shear, d^r,cq/dr, of the equilibrium tide. 
Furthermore, because the daughter energy density is nearly independent of daughter frequency and angular degree, 

'*L i^ ^ weak function of both. 

A. 5. Linear tide coupling coefficient 

In deriving our expression for Kabc, we used the linear homogeneous equations of motion (eqs. |A33| and |A34p to 
simplify the divergence terms in line IA291 Our final expression (lines IA55IIA62)) is therefore only appropriate when 
all three waves are solutions of the homogeneous equations. In order to determine the nonlinear coupling between 
the linear tide ^^^ and a pair of daughters (6, c), we must instead simplify the corresponding divergence terms in line 
IA29I using the inhomogeneous equations of motion; i.e., replace 5(j) ^ 6(j) -\- U in equations (jA33|) and (jA34|) . where 
the normalized tidal potential U{x,t) = — (G'Af/i?)(r/i?)^y£m(0, 0) exp(— iwi). This results in a modified coupling 
coefficient kj^^ = kj^^ -\- kj^J'- the homogeneous part, k^^ , is given by lines IA551IA62I with a replaced by ^jjj^, where 
^ij„ is the solution to equations (JA6IIA8P with U —> U; the inhomogeneous part 



*"= 2MR^ 



^ I drpr' [QThrCr ~ brCh {iAlT + (£ - 1)F} - 5,,c, {3A2t + [l - \)F] 

+bhCh (G + eF)] , (A70) 



where F and G are Fa and Ga with {£a,ma) — >■ {£, m). 

4c and Jbctm (eq. |A23| ). we see that the Ly^a.^,^L^ k,^^^.^ v^l^^.o u^^n ^^-^ s-^ii^v^^-c ^c^l^^^ ,^l v^l^ ^oc^mr^n-^^ 



Comparing kU and Jbdm (eq. |A23| ). we see that the leading order terms bhCh for g-modes cancel in the Jbcim + 2Kl 
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portion of the equilibrium tide coupling coefficient kj^^ (eq. [ID])- Using equation (IA32[) we find 
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(A71) 
(A72) 



where in the second line we used integration by parts and in the third line we made the WKB approximation V • ^ ~ 
^r/^iH (this approximation is well-satisfied throughout the star, including the convective zone). Within the radiative 

zone, the fractional difference between Jtcem and 2k[J is therefore ~ brCr/bhCh ~ [uj/NY ~ 10~^P{^^ and we see 
that there is a large cancellation between these terms. Furthermore, because the b^Cr terms dominate the integrand 

and the radial displacement ^^ is a weak function of period and I (see eq. |A3| ) . Jhctm + 2k\^^ is a weak function of 

period and I. Finally, using iV^ = —g[dlnpldr + l/FiiJ] one can show that the integrand of Jbctm 

in magnitude to k^^ "°'^' (see eq. |A69| ). 



2HilJ is similar 



A. 6. Coefficients of the 'Method 2' amplitude equation 

Here we derive expressions for the coefficients in the method 2 form of the amplitude equation 
We first expand the coefficients in terms of the tidal harmonics: 
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(A75) 



The two other method 2 coefficients, Uab{t) and Kabc, are described in §§ IA.3l and lA.41 respectively. By writing the 
linear displacement as a sum of the equilibrium tide and dynamical tide displacements ^jj„ — $,^^ + ^jyn, we have 
separated coefficients into corresponding components, as indicated by the superscripts "eq" and "dyn" . 

The various coupling coefficients in the definition of Kab are defined in § IA.5I They describe the parametric driving 
of daughter pairs (a, b) considered in §[5l we showed there that the equilibrium tide part of 2Kab cancels strongly with 
Uab and therefore the equilibrium tide part of the sum Uab + 2Kab is much smaller than its individual terms. 

The coefficients Va and Ka describe the nonlinear inhomogeneous driving of a mode a. The term Va corresponds 
to nonlinear tidal driving by the force (^jj,j • V)VC/; one of its components is given by Ja '" — Jab{b — ?> ^dyn) 
where the arrow indicates that the corresponding eigenmode in the expression should be replaced by $,^^^, the 
dynamical tide part of the solution of (jA6IIA8[) with U replaced by the normalized tidal potential U{x,t) — 
— {GM/R){r/RYYim{0,(t))exp{—iujt). The other component Ja is defined similarly but with 1^,^ instead of ^dyn- 
The term Ka describes three- "mode" coupling between mode a and the linear tide and is directly analogous to Kabc 
but with modes b and c at the linear tide values. The components of Ka are given by 

(A76) 
(A77) 
(A78) 

The terms ka and ka ^" ■*" describe driving of mode a by equilibrium tide-equilibrium tide and dynamical 

tide-dynamical tide coupling, respectively, while the term ka~ describes driving by th e cro ss interaction between 
the equilibrium tide and dynamical tide. By the same arguments used to derive equation (|A7ip we have 
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TABLE 1 
COUPLING COEFFICIENTS 



Symbol Coupling of 

Kabc Three eigenmodes H10|l 

ab Equilibrium tide and two eigenmodes 1 140 II 

K* ^'^' Dynamical tide and two eigenmodes I I41II 

^jeq eqj Equilbrium tide to itself and an eigenmode l)A81[l 

i^a '*'" ''^ Dynamical tide to itself and an eigenmode IIA82I I 

^i^eq yn; Equil. tide, dyn. tide, and an eigenmode IIA83II 



Note. — Numbers in parentheses refer to the equation where 
the coefficient is defined. All k's not equal to K^hc a-re simply 
versions of K^hc with c and/or b equal to components of the linear 
tide. 

The equivalent term involving the equilibrium tide can be simplified further because V • ^^„ — 
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91nr 
drr'+^5palr,.n- (A80) 



As in equation (jA7ip . we find that there is significant cancellation between J a and 2k ^^^{b — >■ ^^y^) and like- 
wise for the equivalent equilibrium tide expression. It is therefore convenient to define the three types of nonlinear 
inhomogeneous driving coefficients as follows: 

^^cq-oq) ^ j{cq) ^ ^(^cq-eq) ^^gj^j 

^(^dyn-dyn) ^ j^dyn) _^ ~^dyn-dyn) ^^g2) 

^(^cq-dyn)^2Ki'=<l-dy"). (A83) 

The properties of these three coefficients are described in § IB.3I 

B. PROPERTIES OF THE NONLINEAR COUPLING COEFFICIENTS 

In this appendix we describe the properties of the nonlinear coupling coefficients for a solar model. For reference, 
we list the various coupling coefficients in Table [H We first discuss the two coupling coefficients corresponding to 

parametric driving: the dynamical tide coefficient k\^/'^ (iJ IB.ip and the equilibrium tide coefficient kj^J^' (^ \B.2\i . We 

show that nonlinear driving by the dynamical tide Tf^/"' oc kj^J^^ and equilibrium tide Fj^'^ ex kj^^ are, in general, 
highly localized, with most of the coupling occurring at a specific location in the star: near the core in the case of 

^bc'^ and at the radiative-convective interface in the case of Fj^"'^ . The only exception is equilibrium tide driving of 
self-coupled daughters, which tends to be more global and, as a result, can potentially lead to driving throughout the 

star. 

In ? IB. 31 we describe the coupling coefficients corresponding to nonlinear inhomogeneous driving Va and Ka- Just 
as the linear inhomogeneous term Ua can be expressed in terms of a dimensionless coefficient laim times a tidal factor 
~ e (see eq. |14]'). the nonlinear inhomogeneous term Va + Ka can be expressed in terms of dimensionless coupling 
coefficients times a tidal factor ^ e^ . By expressing the linear tide as the sum of the uj — equilibrium tide and 
the dynamical tide, ^u^ = ^g^ -I- ^^ynj these coefficients can be broken into three types of interactions: equilibrium 

tide-equilibrium tide coupling Ka °'^ , dynamical tide-dynamical tide coupling Kq ■*'" ^'^ , and the cross coupling 

between the equilibrium tide and dynamical tide Ka~ ■ We define these coefficients in ij IA.6I (egs. |A81| . |A82| . 

and |A83j ). In general, we find that for a solar binary eKa~ is much smaller than laim- However, both e^i ^"~ ^'^' 

and EKa ■*" can be much larger than laim for certain modes (including p-modes). As we discuss in §[51 this implies 
that these modes are nonlinearly driven to energies that are much larger than their linear values. 

In order to calculate these various coefficients we need to know the linear tidal response, as described by the 
equilibrium and dynamical tides ^^q and ^dyu- We describe how we calculate these in ij lA.ll 



B.l. Dynamical tide coupling coefficient kj^J''^ 

PP 
mode coupling coefficient Kabc when mode a is the parent eigenmode most linearly resonant with the tide. This is 



The spatial structure of the dynamical tide coupling coefficient kj^J"^ is approximately equal to that of the three 
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Fig. 12. — Left panel: Radial profile of the tliree-mode coupling coefficient k^^i,,. for a solar model. The parent mode is {ia,Tna,na) = 
(2, 0, 486) and the daughter is a self-coupled mode {b = c) with (^j,, m;,, n;,) = (2, 0, 971). The periods of these modes are Pa ~ ^i,/2 ~ 5 days. 
Right panel: |Ka6cl ^^ ^ function of An = ni, — ric- The parent mode a and daughter mode b are the same as in the left panel. The other 
daughter mode, whose radial order rtc we vary, has {£c, "mc) = (2, 0). 

because the dynamical tide is typically dominated by this single mode (e.g., the mode with n = 486 in Fig. [T]). The 

approximate magnitude of kU'^ can be found by multiplying Kate by 0JaIae.m/{^a — ila) (see eq. [21])- Since it is 
more common to express three- wave coupling coefficients in terms of three eigenmodes, we show results for Kabc rather 

than i%c^ ■ 

In the left panel of Figure [12] we show the integrand of Kabc for a linearly resonant parent coupled to a self-coupled 
daughter b = c. The coupling coefficient Kabc involves an integral over an odd number (three) of high-order modes. 
The integrand peaks in the core since the radial and horizontal displacements i^r,^h) of high order modes increase 
with decreasing r. For a standing wave, nearly all the contribution to Kabc comes from a small region between the 
daughter's inner turning point where ujb — N and the parent's inner turning point where uja — N] since ujb < uJa 
and dN/dr > 0, the parent's turning point lies above the daughter's. In this region, the parent is evanescent and 
the displacements £,a,r,£,a,h oc r^"~^ do not oscillate. As long as the factor cos [(fcb — kc)r] is roughly constant in this 
region, the modes add coherently and can exchange energy through nonlinear interactions. 

In the right panel of Figure [1^ we show \Kabc\ as a function of An = Ub — Uc for £ — 2 modes. The radial order 
of the parent mode a is Ua = 486. For |An| < Ua, we find that \Kabc\ is approximately constant and close to its 
maximum value (to within a factor of a few); otherwise it is much smaller than the maximum value. This is because 
in the propagation zone for each wave, the integrand is proportional to a factor cos(fcaf) cos(fcbr) cos(fcc''), where ka 
is the radial wavenumber of mode a, etc. This product of three cosines can be combined into a sum of terms of the 
form cos [{ka i kb ± kc) r]. The coupling coefficient is thus maximized when fc^ ± fc;, ± fee ~ 0, the us ual condition for 
momentum conservation. For high-order modes this is roughly equivalent to jn^ — nd ^ria (see also lWu fc GoldreichI 
1200 It) . Hence although there is no rigorous selection rule for the radial direction, as there is for £ and m, there is an 
approximate selection rule. When momentum conservation is satisfied in the appropriate region, the magnitude of 
Kabc is near its maximum value. When it is not satisfied, the coupling coefficient is much smaller than this maximum 
possible value. 

The approximate selection rule |nf, — nd ^ Ua implies that there can be a very large number daughter pairs with 
similar values of Kabc- As a result, each daughter couples well to not only one other daughter but to N ^ Ua other 
daughters. As we show in § 17.21 this implies that collective driving can be very important for the dynamical tide. 

B.2. Equilibrium tide coupling coefficient Kj^^ 

We show the radial dependence of k^^' for the equilibrium tide coupled to a self-coupled daughter mode & = c in 
Figure [13] and for daughter modes b ^ c in Figure [14] Unlike Kabc, the integrand of the equilibrium tide coupling 

coefficient kj^J^ contains a product of an even number of modes. More precisely, since the displacement due to the 
equilibrium tide increases monotonically with increasing radius and there is a factor of cos [{kb ± kc)r] in the integrand, 
momentum conservation requires that locally \kb\ ~ \kc\. The coupling is therefore maximized for self-coupled modes 
6 = c; as the number of nodes differs by | An| — \nb ~ nd =0,1,2..., the equilibrium tide coupling decreases rapidly 
due to cancellation between the modes. As a result, unlike with the dynamical tide, collective driving does not appear 
to be important for the equilibrium tide. 



NONLINEAR TIDES IN CLOSE BINARY SYSTEMS 



35 



0.02 



A 



0.01 










r / R 



0.05 



0.04 



^^ 0.03 



0.02 - 



0.01 - 




t=2 



10 



15 20 



Pb (days) 



,(cq) 



Fig. 13. — Left panel: Radial profile of the equilibrium tide coupling coefliicient kJ,^ for a solar model. The parameters of the equilibrium 
tide are {£,m) = (2,0) and the daughter is a self-coupled mode (6 = c) with (^6,ni6,»i6) = (2,0,971), which corresponds to a period of 
P(, ~ 10 days. The individual oscillations in the integrand dj/dlnr are too closely spaced to be seen well at the resolution of the figure. 



Right panel: Kj '^ as a function of period P(, for self-coupled daughters b = c with m = mi, = and 
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Fig. 14. — Left panel: Radial profile of the equilibrium tide coupling coefficient f^l^^ for the {£, m) = (2, 0) equilibrium tide coupled to 
different daughter pairs (fe, c). The four curves are for {£i,,rni,,ni,) = (2,0,971) and {Ic'm-cnc) = (2,0, [971,972,973,975]). Right panel: 



Hj^^ as a function of An = U), — Uc for ii, = £c = "2, "m = m,;, 



We also find that for self-coupling, k/'^' 



: 0, and Pfc 



10 days. 



T,j^ is not strongly concentrated in radius. However, as |An| increases the 
coupling becomes increasingly concentrated in the region just below the radiative-convective interface. As shown in 

the right panel of Figure [131 i^bc ^^ ^^^° ^ weak function of daughter period and £. 



B.3. 



Inhomogeneous coupling coejficients 
In Figure[T2]we show the magnitudes of the three nonlinear inhomogeneous coupling coefficients Ka , 

(cq-dyn) 



^(dyn-dyn) 



and Ka ■ These represent driving of a mode a due to equilibrium tide-equilibrium tide coupling, dynamical tide- 

dynamical tide coupling, and equilibrium tide-dynamical tide coupling, respectively. The modes in the left panel range 
from high-order g-modes with n sa 10'^ to low-order p-modes with n < 30. Since the acoustic cutoff of the solar 
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Fig. 15. — The nonlinear inhomogeneous coupling coefficients kJ, , ta , and kJj as a function of mode frequency uia 

in units of loq = {GM/ B?')^''^ for a solar model. The daughters are all (^a,m,a) = (2,0) modes coupled to the {i,m) = (2,0) dynamical 
tide at a P = 10 day driving period. Those in the left panel range from high-order g-modes to low order p-modes while those in the right 



panel are high-order p-modes (700 < ria < 1400). In the right panel we divide kJ, 

(oq-dyn) 



(dyn — dyn) 



by 10 in order to show it on the same scale as 



atmosphere is ss BOwq, these p-modes are trapped. The right panel shows results for high-order p-modes that are all 

well above the acoustic cutoff and therefore do not reflect at the solar surface. 

The integrand of the coupling coefficient Ka is similar to that of the linear overlap laim (compare eqs. |A15| 



(oq-oq) 



elai 



and |A80| ). They therefore have very similar dependences on Wq (cf. Fig. [TlT) and for £ = 2 we find Ka 

Driving by Ha , which 

compared to linear driving 

The coupling coefficient Ka ■^" '*'" is by far the largest of the three coefficients, typically by at least 2-3 orders of 
magnitude. From scaling arguments we see that its magnitude is roughly equal to Kabb{Ib^b/ ^b)^ , where mode b here 

refers to the dynamical tide mode and Kabb ~ 5 x IQ^P^^^ (eq. [^67]). We find that Ka '^"^ "*'" is large as long as 
the frequency of mode a (the mode being driven) is larger than the frequency of the dynamical tide. This is because 
there is then a region below mode a's inner turning point where a is evanescent and the integrand does not oscillate 



Driving by Ka , which draws all of its energy from just the equilibrium tide, therefore appears to be insignificant 



about zero. This behavior is identical to that of 



,(dyn) 
"■be 



('iJ IB.ip except that there the dynamical tide serves as the high 
frequency mode in the triplet whereas here mode a serves as the high frequency mode. For the case of a solar binary, 

£Ka ■*'" ■*'" is larger than lahn for g-modes with periods Po ^ 1 day- This suggests that inhomogeneous driving by 
the dynamical tide may be important for these modes. We discuss this possibility in § [HI 

We find that Ka '^" ^" is especially large for high-order p-modes and increases with increasing uja (right panel of 
Figure [T5|) . Physically, this is because high-order p-modes penetrate deeply into the solar core, where the amplitude 
of the dynamical tide peaks. The p-modes turning point occurs where LOa equals the Lamb frequency Si^{r) = AaCg/r. 
The sound speed is nearly constant in the core (cs ~ 5 x 10^ cm s~^) and high-order ^^ = 2 p-modes do not reflect 
until they reach a depth r/R w 0.003 (wa/ 10'^ wq)^^. In §|n]we briefly discuss whether these p-modes, which cannot 
form standing waves since they are above the acoustic cutoff, can be driven to significant amplitudes. 

The coupling coefficient Ka ^^ is significantly larger than Ka "'^ for high-order g- and p-modes. For high-order 
g-modes, both coefficients are dominated by coupling within the convection zone as all three waves in the triplet are 
long wavelength in the convection zone and can thus couple well there. However, unlike Ka , we find that Ka '"^ 

is weakly dependent on uja for g-modes. This is because the integrand of the latter depends on {^r,a,^h,a), both of 

which vary weakly with uja within the convection zone, where the mode is evanescent. By contrast, Ka depends 

on 6 Pa (see eq. |A80p which is a strong function of uJa within the convection zone. 

For high-order p-modes, the coupling Ka '^' does not occur within the convection zone since the p-modes are 
oscillatory there. Instead, the coupling occurs in the radiative interior where the p-mode wavelength Ap « 2'KCs/uJa is 
comparable to the dynamical tide wavelength Adyn ~ 2'KUjr / KN . Since Ap decreases monotonically with r and Adyn 
increases monotonically with r, there is a region where the two wavelenths cross for p-modes with uJa in a certain 
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Fig. 16. — Radial profile of Ka for the high order p-mode (la = 2, ria = 1200, nia = 0) coupled to the (£, m) = (2, 0) equilibrium 

tide and the (£, m) = (2, 0) dynamical tide at a P = 10 day driving period. The dashed lines show the variation of the radial wavelength 
of the p-mode Ap and the dynamical tide A;jy„ in units of R. 



range. In Figure [TB] we show the radial profile of Ka 



,(cq-dyn) 



over such a region for an Ua = 1200 p-mode. The integrand 



of Ka '^' oscillates symmetrically about zero in regions where the wavelengths are very different. However, at 
r ~ 0.37R the two wavelengths cross (dashed lines) and the two modes couple coherently over a region of size sa 0.05i?. 
We find that over the frequency range lO^Ps < Wq/wq '^ W^P^ there is a radius in the radiative interior of a solar-type 
star where Adyn = Ap ; a small fraction of this frequency range is shown in the right panel of Figure [151 

C. INTERACTION ENERGY AND ORBITAL EVOLUTION 

Here we derive the acceleration of the secondary due to the oscillation modes excited in the primary assuming the 
modes are all standing waves. This is most easily accomplished by plugging the tidal potential (eq. |11| ) into the 
interaction Hamiltonian (eq. [2]), and using the integrals defined in equations p2)) and (fT3| . We find 



Hint — ^ 



^ GMM'WimR'^ 
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Z)i+i 



fa + l^Yl '^'^' 



Mmqalb 



In the center of mass frame, this interaction energy gives the following radial and angular accelerations: 



1 dHint 
fi oD 

1 9i?int 



G{M + M') 



aim 



a$ : 



fiD a$ 



G{M + M') ^ . fR 

aim 




-im^ 



(Cl) 



(C2) 



(C3) 



where ^ — MM' /{M + M'). The expressions in equations (IC2[) and (|C3|) may also be derived from the perturbed 
potential of the prim ary evaluated at the position of the secondary. In doing so, the nonlinear expression for 5p from 
iSchenk et al.l ()2002[ ) must be used to compute the quadrupole moment. 

For orbital frequen cy f? = \G(M +M') / a?]^/'^ , the radial and angular accelerations due to the tide imply the following 
changes in the orbit (jMurrav fc Dermottll2000l ): 






[aoe sin $ -I- a$ (1 -I- e cos $)] 



Q.G 



[ao sin $ + a$ (cos $ + cos E)] 



(C4) 
(C5) 



where cosE = (e -I- cos$)/(l + ecos$). Writing the mode amplitude in amplitude-phase form g^ = AaC ""^"^ , where 
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Aa and (pa are real, the decay and circularization rates (eqs. |C4| and \C5\ ) can be written 



d/a 
e/e 



SI J2 WemX, 



em 

k 



aimk 



h 



lalmAa sin ((^a - fcfii) 



" 9 X] "^abfm^aA sin {kpa +(/?& + fef^i) 



(C6) 



where /i = 2fc and /2 = (e ^ — l) [/c — ni/Vl — e^J . At resonances the arguments in the sines can be time independent 
with a constant lag angle that depends on the damping rates and detunings (see, e.g., eq. [25] and Appendix [P]) . 

D. NONLINEAR EQUILIBRIUM FOR A SIMPLE THREE MODE SYSTEM 

In this appendix we solve for the nonlinear equilibrium (i.e., saturation) of a simple 3 mode network in which a single 
parent a, driven by the linear tide, is coupled to a pair daughters {b,c). We assume that the daughters are distinct, 
not self-coupled, and not driven by the linear or nonlinear tide. The amplitude equations for this system are thus 



qa + (iuJa + la)qa = i^a [Ua{t) + ^Kqlq*] , 

qb + {iuJb + "tb)qb = '2i(^bKql<lcy 
qc + {iujc + jc)qc = '2iuJcK.qlql, 



(Dl) 
(D2) 
(D3) 



where k = Habc and the factors of two come from assuming that the modes are distinct from one another. Focusing on 
a particular harmonic uj such that Ua{t) = t/ae~"^*, let qa = QQ,e~'^'^°~^°-'*, where a = {a,b,c} and Qa is a constant 
complex amplitude. If Aq = iOa — uj and Ab + Ac = lu + LOb + LUc the time dependences drop out yielding the equilibrium 
solution 



Aae'^-^ {iAa+-fa)=i^a 



Ua + 2KAbAce 



~i(5i,+'5c) 



Abe'"' {lAb 



lb. 



= 2iiUbKAaA,e-'^^-+^^\ 



A^e^'^ (iAc + jc)=2iu:cKAaAbe-'^'^+''\ 



(D4) 

(D5) 
(D6) 



where we have written Qa — Aac'^" with Aa and 6a realo 

To solve for the parent's nonlinear equilibrium energy, divide the two daughter equations (IDSP and (|D6p to get 






A^ _ Ub f AbAc + 7;, 7c - ^(Afc7c - Ac7fc) 



(D7) 



Wc V A^ + % 

Since the left hand side is real, the imaginary part on the right hand side must vanish, Af,/Ac — Jb/lc^ and we find 



^b 



Al 



UJclb 



Multiplying the two daughter equations we obtain 



Ai 



Iblc 



AK-^UJbiOi 



[(1 -^2)cos2(5-2/isin2(5 + i(2^cos2(5+ (1 - /i^) sin 2.5)] , 



(D8) 



(D9) 



where ^ = (Af, + Ac)/(76 + 7c) and 6 = Sa + 6b + 6c. Since the left hand side is real, the imaginary part must vanish 
so that 

tan 25= ^. (DIO) 

1 — /x^ 

There is a sign ambiguity which is resolved by noting that A^ > 0; if, for example, < fj, < 1, choose the quadrant 
such that sin 26 = 2/x/(l + /i-^) > and cos 26 — —(1 — M'^)/(l + M^) < 0- We thus find the parent's equilibrium energy 

Iblc 






An'^uJbUJi 



(1 + M^) 



(Dll) 



which we see equals its threshold energy Eth (eq. [SD]). 

To solve for the daughters' nonlinear equilbrium energy, rearrange the parent's equation (|D4p and use equations 
(|D8l) and ((DTT|) to get 



^b ■ 



7c\/TT7? 



Aa COS 6 + 7a sin (5 — ^a—f- cos{6 — 6a) + i I Aa sin (5 — 7a cos (5 — Wa-p sm{6 — 6a) 

Aa \ Aa 



An'^UJaU!, 

lerived thf 
re |r| = i 
unstable if the parent's linear energy E^ Un > ^thi where Eth is 



(D12) 



I n § \E\ we derived the stability criteria for this type of system 
(eq. 1371 V Here |r| = 2KAay/uJi,cJc and the linear tidal flow is 



given by equation ((50). 
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Since A^ is real, the imaginary part on the right hand side must vanish 



WQ-pSm(d- 



6n 



Aa sin S -"fa COS S, 



(D13) 



which, with equations (jD10[) and (IDllI) and the definition oi S, determines (5a, (5;,, and Sc- The daughters' equihbrium 
solution is thus given by 
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- 4^ - 
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AK-^UJai^c 



AaC0S(5 + 7a sin (5 - uja-f- cos{6 ~ Sa) 



(D14) 



and Ec,cq 
{seem, 



Eb.cq^db/'-^bJc- In the limit -Ea,iin 3> -Ea,oq, as is appropriate for solar binaries and hot Jupiter systems 



Eb,cq — 



Jci^b 



Ua 

2k 



(D15) 



7fcWc 
In the limit i?a, lin -^ Eth — Ea^cq the term in brackets in (|D14p vanishes and Ab,Ac — > 0. 

E. LOCAL NONLINEAR GROWTH RATE 

In this Appendix we calculate the local nonlinear growth rate due to three mode coupling and nonlinear tidal driving. 
The local rate per unit volume at which the nonlinear interactions do work is given by the fluid velocity dotted into 
the nonlinear force: e{x) — 2^ • (/2[^!^] — p(^ • V)VC/). For a self-coupled daughter wave b 



eb{x)^2\qb\^, ■ {\qa\\qb\f2[^a,^b] " PMi^b ' V)VL/) 



^m\ 



19a I 



Jbb 



drdQ drdfl 



Eo, 



(El) 



where we assume that the interactions drive the daughter at half the tidal frequency lj/2. In the WKB approximation, 
the local energy density of an internal gravity wave is 



91- 



eW = -jT 



d^E 
drdfl 



dn^pN^lqfCEo. 



(E2) 



The local nonlinear growth rate of the daughter over a region Ar much smaller than r (and the scale over which the 
background quantities vary) but much larger than its radial wavelength 2Tr/kr^b is therefore 



r.(r): 
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Ar 
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eb{r) 

Eq 
dE/dr 



dr 



, dKr 
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bb 



dr 



dUbb 
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(E3) 



The product of the phases of the traveling waves (e.g., exp[i(0a + 2(j)b)]) is constant by momentum conservation and 
we have integrated over angles since the horizontal group velocity is faster than the vertical group velocity by a factor 
of kr/kh = N/u j > 1 . 
By equation (|A65p . the local growth rate of daughters due to three- wave coupling with a parent a is thus 



b.K - 



^\qa\ 



r+Ar 



Eg dK, 



abb 



dE/dr dr 



dr ~ ujT\qkrCr\a - ^T\qAkhih\a- 



(E4) 



We thus see that for the coupling of high-order g-modes, three-wave interactions act as a Kelvin-Helmholtz-like shear 
instability that is driven by the parent's horizontal motion ^^; the growth rate is given by the oscillation frequency times 
the strength of the horizontal shear qAkh£,h- Whereas overturning of the vertical stratification requires If^l/cr^i. > 1, a 
g-mode is unstable to three-wave coupling even for \q\kr£,r ^ 1, although at small enough amplitude the growth rat e 
of the daughters becomes smaller than their linear damping rate 76 (see, e.g.. lDrazinlll977tlSonmor fc Klaassenlll997l ). 
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